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PREFACE. 



Theory of relativity is a physical discipline which arose in 
the begiiming of XX-th century. It has dramatically changed 
traditional notion about the structure of the Universe. Effects 
predicted by this theory becomes essential only when we describe 
processes at high velocities close to light velocity 

c = 2.998 • Iff'km/sec. 

In XlX-th century there was the only theory dealing with such 
processes, this was theory of electromagnetism. Development of 
theory of electromagnetism in XlX-th century became a premise 
for arising theory of relativity. 

In this book I follow historical sequence of events. In Chapter I 
electrostatics and magnetostatics are explained starting with first 
experiments on interaction of charges and currents. Chapter II is 
devoted to classical electrodynamics based on Maxwell equations. 

In the beginning of Chapter III Lorentz transformations are 
derived as transformations keeping form of Maxwell equations. 
Physical interpretation of such transformation requires uniting 
space and time into one four-dimensional continuum (Minkowsky 
space) where there is no fixed direction for time axis. Upon 
introducing four-dimensional space-time in Chapter III classical 
electrodynamics is rederived in the form invariant with respect to 
Lorentz transformations. 

In Chapter IV variational approach to describing electromag- 
netic field and other material fields in special relativity is con- 
sidered. Use of curvilinear coordinates in Minkowsky space and 
appropriate differential-geometric methods prepares background 
for passing to general relativity. 
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In Chapter V Einstein's theory of gravitation (general rela- 
tivity) is considered, this theory interprets gravitational field as 
curvature of space-time itself. 

This book is addressed to Math, students. Therefore I paid 
much attention to logical consistence of given material. References 
to physical intuition are minimized: in those places, where I 
need additional assumptions which do not follow from previous 
material, detailed comment is given. 

I hope that assiduous and interested reader with sufficient 
preliminary background could follow all mathematical calculations 
and, upon reading this book, would get pleasure of understanding 
how harmonic is the nature of things. 

I am grateful to N. T. Ahtyamov, D. I. Borisov, Yu. P. Ma- 
shentseva, and A. I. Utarbaev for reading and correcting Russian 
version of book. 



November, 1997; 
November, 2003. 



R. A. Sharipov. 



CHAPTER I 



ELECTROSTATICS AND MAGNETOSTATICS 



§ 1. Basic experimental facts and unit systems. 



Quantitative description of any physical phenomenon requires 
measurements. In mechanics we have three basic quantities and 
three basic units of measure: for mass, for length, and for time. 



Quantity 


Unit 


Unit 


Relation 




in SI 


in SGS 


of units 


mass 


kg 


9 


1^5 = 10^5 


length 


m 


cm 


lm = 10^ cm 


time 


sec 


sec 


1 sec = 1 sec 



Units of measure for other quantities are derived from the 
above basic units. Thus, for instance, for measure unit of force 
due to Newton's second law we get: 

(1) N ^ kg ■ m - sec'"^ in SI, 

(2) dyn = g ■ cm ■ sec~'^ in SGS. 

Unit systems SI and SGS are two most popular unit systems 
in physics. Units for measuring mechanical quantities (velocity, 
acceleration, force, energy, power) in both systems are defined in 
quite similar way. Proportions relating units for these quantities 
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can be derived from proportions for basic quantities (see table 
above). However, in choosing units for electric and magnetic 
quantities these systems differ essentially. 

Choice of measure unit for electric charge in SGS is based on 
Coulomb law describing interaction of two charged point. 



Coulomb law. Two charged points with charges of the same 
sign are repulsing, while points with charges of opposite signs are 

attracting with force proportional to 
quantities of their charges and in- 
verse proportional to square of dis- 
tance between them: 



■o 



■0 ©■ 



(1.1) Fr^^ 



► Measure unit of charge in SGS is 

chosen such that coefficient in for- 
mula (1.1) is equal to unity. Hence 
^'''9- we have the following relation: 

unit of charge in SGS = dyn^^^ ■ cm= g^l^ ■ cn^^^ ■ sec~^. 



Coulomb law itself then is written in form of the equality 




Force F defined by the relationship (1.2) is very strong. How- 
ever, in everyday life it does not reveal itself. This is due to the 
screening. The numbers of positive and negative charges in nature 
are exactly balanced. Atoms and molecules, which constitute all 
observable matter around us, have the same amount of positive 
and negative charges. Therefore they are electrically neutral in 
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whole. Force (1.2) reveals itself in form of chemical links only 
when atoms are pulled together. 

Electric current arises as a result of motion of charged points. 
This occurs in metallic conductor, which usually have lengthy 
form (form of wire). Current in such conductor is determined by 
the amount of charge passing through it within the unit of time. 
Therefore for unit of current we have: 

unit of current in SGS = unit of charge in SGS ■ sec~^ = 
= g'/^ ■ cm^/^ ■ sec-\ 

Let's consider straight conducting rod of the length I. Current 
in it leads to misbalance of charges in its ends. Charges of definite 
sign move to one end of the rod, while lack of these charges in 
the other end of the rod is detected as the charge of opposite 
sign. Then Coulomb force (1.2) arises that tends to recover 
balance of charges in electrically neutral rod. This means that 
in such rod current could not flow in constant direction during 
long time. Another situation we have with conductor of the 
form of ring or circuit. Here current does not break the balance 
of charges. Direct current can flow in it during unlimitcdly 
long time. Circular conductor itself thereby remains electrically 
neutral and no Coulomb forces arise. 

In spite of absence of Coulomb forces, in experiments the inter- 
action of two circular conductors with currents was detected. This 
interaction has other nature, it is not due to electrical, but due to 
magnetic forces. The magnitude of magnetic forces depends es- 
sentially on the shape and mutual arrangement of circular conduc- 
tors. In order to reveal quantitative characteristics for magnetic 
forces one should maximally simplify the geometry of conduc- 
tors. For this purpose they are deformed so that each possesses 
straight rod-shaped part of suflaciently big length I. These rod- 
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shaped parts are arranged par- 
allel to each other with the dis- 
tance r between them. In the 
limit, when I is much larger than 
r, this configuration of conduc- 
tors can be treated as a pair 
of infinitely long parallel conduc- 
tors. In experiments it was found 
that such conductors do interact 
according to the following law. 

Ampere law. Force of inter- 
action of two infinite parallel 
conductors with currents per unit 
length of them is proportional ti the values of currents in them 
and inverse proportional to the distance between them: 





Fig. 1.2 



(1.3) 



F 
1 



hh 



Two co-directed currents attract each other, while opposite di- 
rected currents repulse each other. 

The \init of current in SGS was already introduced above. 
Therefore coefficient of proportionality in formula (1.3) is unique 
quantity that should be determined in experiment. Here is the 



inverse square of velocity. Therefore formula (1.3) in SGS is 
written as 



(1.4) 



F 
1 



2 hh 



Constant c in (1.4) is a velocity constant. The value of this 
constant is determined experimentally: 

(1.5) c w 2.998 • 10^° cm/ sec. 
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As we shall see below, constant c in (1.5) coincides with velocity 
of light in vacuum. Numeric coefficient 2 in (1.4) is introduced 
intentionally for to provide such coincidence. 

In SI measure unit of current 1 A (one ampere) is a basic unit. 
It is determined such that formula (1.3) is written as 

Here tt = 3.14. . . is exact (though it is irrational) mathematical 
constant with no measure unit. Constant /jq is called magnetic 
susceptibility of vacuum. It has the measure unit: 

(1.7) = Att ■ IQ-^ N ■ A-"^ . 

But, in contrast to constant c in (1.5), it is exact constant. Its 
value should not be determined experimentally. One could choose 
it to be equal to unity, but the above value (1.7) for this constant 
was chosen by convention when SI system was established. Due 
to this value of constant (1.7) current of 1 ampere appears to 
be in that range of currents, that really appc^ar in indiistrial and 
household devices. Coefficient An in denominator (1.6) is used in 
order to simplify some other formulas, which arc more often used 
for engineering calculations in electric technology. 

Being basic unit in SI, unit of current ampere is used for 
defining unit of charge of 1 coulomb: IC = 1^ • Isec. Then 
coefficient of proportionality in Coulomb law (1.1) appears to be 
not equal to unity. In SI Coulomb law is written as 

(1.8) i.= J_^. 
^ ' 47reo 

Constant eo is called dielectric permittivity of vacuum. In con- 
trast to constant /xq in (1-7) this is physical constant determined 
experimentally: 

(1.9) eo » 8.85 • IQ-^^ • TV"! • m'^. 
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Constants (1.5), (1.7), and (1.9) are related to each other by the 
following equality: 

(1.10) c = fti 2.998 • 10^ ml sec. 

V eo Mo 

From the above consideration wc see that SGS and SI systems 
differ from each other not only in the scale of units, but in for- 
mulas for two fundamental laws: Coulomb law and Ampere law. 
SI better suits for engineering calculations. However, derivation 
of many formulas in this system appears more huge than in SGS. 
Therefore below in this book we use SGS system. 

Comparing Coulomb law and Ampere law we sec that electrical 
and magnetic forces reveal themselves in quite different way. 
However, they have common origin: they both are due to electric 
charges. Below we shall see that their relation is much more 
close. Therefore theories of electricity and magnetism are usually 
united into one theory of electromagnetic phenomena. Theory of 
electromagnetism is a theory with one measurable constant: this 
is light velocity c. Classical mechanics (without Newton's theory 
of gravitation) has no measurable constants. Newton's theory of 
gravitation has one constant: 

(1.11) 7 K. 6.67 • 10"^ cm^ • 5"^ • sec''^. 

This theory is based on Newton's fourth law formulated as follows. 

Universal law of gravitation. Two point masses attract 
each other with the force proportional to their masses and in- 
verse proportional to the square of distance between them. 

Universal law of gravitation is given by the same formula 



(1.12) 



_ M1M2 
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in both systems: in SGS and in SI. 

According to modern notion of nature classical mechanics and 
Newton's theory of gravitation are approximate theories. Cur- 
rently tlicy arc replaced by special theory of relativity and general 
theory of relativity. Historically they appeared as a result of 
development of the theory of electromagnetism. Below we keep 
this historical sequence in explaining all three theories. 

Exercise 1.1. On the base of above facts find quantitative re- 
lation of measure units for charge and current in SGS and SI. 

§ 2. Concept of near action. 

Let's consider pair of charged bodies, which are initially fixed, 
and let's do the following mental experiment with them. When 
we start moving second body apart from first one, the distance r 
begins increasing and consequently force of Coulomb interaction 
(1.2) will decrease. In this situation we have natural question: 
how soon after second body starts moving second body will feel 
change of Coulomb force of interaction? There are two possible 
answers to this question: 

(1) immediately; 

(2) with some delay depending on the distance between bodies. 

First answer is known as concept of distant action. Taking this 
concept we should take formula (1.2) as absolutely exact formula 
applicable for charges at rest and for moving charges as well. 

Second answer is based on the concept of near action. Ac- 
cording to this concept, each interaction (and electric interaction 
among others) can be transmitted immediately only to the point 
of space infinitesimally close to initial one. Transmission of any 
action to finite distance should be considered as a process of 
successive transmission from point to point. This process always 
leads to some finite velocity of transmission for any action. In 
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the framework of the concept of near action Coulomb law (1.2) is 
treated as approximate law, which is exact only for the charges 
at rest that stayed at rest during sufficiently long time so that 
process of transmission of electric interaction has been terminated. 

Theory of electromagnetism has measurable constant c (light 
velocity (1.5)), which is first pretender for the role of transmission 
velocity of electric and magnetic interactions. For this reason 
electromagnetic theory is much more favorable as compared to 
Newton's theory of gravitation. 

The value of light velocity is a very large quantity. If wc settle 
an experiment of measuring Coulomb force at the distances of the 
order of r « 10 cm, for the time of transmission of interaction we 
would get times of the order of t « 3 • 10^^" sec. Experimental 
technique of XlX-th century was unable to detect such a short 
interval of time. Therefore the problem of choosing concept could 
not be solved experimentally. In XlX-th century it was subject for 
contests. The only argument against the concept of distant action 
that time, quite likely, was its straightness, its self-completeness, 
and hence its scarcity. 

In present time concept of near action is commonly accepted. 
Now we have the opportunity for testing it experimentally in the 
scope of electromagnetic phenomena. Let's study this concept 
more attentively. According to the concept of near action, process 
of transmitting interaction to far distance exhibits an inertia. 
Starting at one point, where moving charge is placed, for some 
time this process exist in hidden form with no influence to both 
charges. In order to describe this stage of process we need to 
introduce new concept. This concept is a field. 

Field is a material entity able to fill the whole space and able 
to act upon other material bodies transmitting mutual interaction 
of them. 

The number of fields definitely known to scientists is not big. 
There are only four fundamental fields: strong field, weak field, 
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electromagnetic field, and gravitational field. Strong and weak 
fields are very short distance fields, they reveal themselves only 
in atomic nuclei, in collisions and decay of elementary particles, 
and in stellar objects of extremely high density, which are called 
neutron stars. Strong and weak interactions and fields are not 
considered in this book. 

There arc various terms using the word field: vector field, tensor 
field, spinor field, gauge field, and others. These are mathematical 
terms reflecting some definite properties of real physical fields. 



Let's apply concept of near action to Coulomb law for two 
charged points. Coulomb force in the framework of this concept 
can be interpreted as follows: first charge produces electric field 
around itself, and this field acts upon other charge. Result of such 
action is detected as a force F applied to second charge. Force 
is vectorial quantity. Let's denote by F vector of force and take 
into account the direction of this vector determined by verbal 
statement of Coulomb law above. This yields 



Here ri and r2 are radius-vectors of points, where charges Qi and 
Q2 are placed. Let's consider vector E determined as the ratio 
E = F/Q2- For this vector from formula (3.1) we derive 



§ 3. Superposition principle. 



(3.1) 



F = Qi Q2 



r2 - ri 



r2 - ri 



(3.2) 



r2 - ri 



|r2 - ri|3" 



Vector E depends upon the position of first charge and upon 
its value. It depends also on the position of second charge, but 
it doesn't depend on the value of second charge. One can take 
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vector E for quantitative measure of electric field produced by first 
charge Qi at the point r2, where second charge is placed. Vector 
E can be determined by formula (3.2) or it can be measured 
experimentally. For this purpose one should place test charge q 
to the point r2 and one should measure Coulomb force F acting 
upon this test charge. Then vector E is determined by division of 
F by the value of test charge q: 



Now consider more complicated situation. Suppose that 
charges Qi, . . . , Qn are placed at the points ri, . . . , r„. They 
produce electric field around them, and this field acts upon test 
charge q placed at the point r. This action reveals as a force F 
applied to the charge q. Again we can define vector E of the form 
(3.3) and take it for the quantitative measure of electric field at 
the point r. This vector is called vector of intensity of electric 
field or simply vector of electric field at that point. 

Generally speaking, in this case one cannot be a priori sure 
that vector E does not depend on the quantity of test charge q. 
However, there is the following experimental fact. 

Superposition principle. Electric field E at the point r 
produced by a system of point charges Qi, . ■ . ,Qn is a vectorial 
sum of electric fields that would be produced at this point by each 

charge Qi, . . . , Qn separately. 

Superposition principle combined with Coulomb law leads to 
the following formula for the intensity of electric field produced 
by a system of point charges at the point r: 



(3.3) 



E = F/q. 



n 



(3.4) 




§3. SUPERPOSITION PRINCIPLE. 



17 



Using superposition principle, one can pass from point charges 
to continuously distributed charges. Suppose that the number of 
point charges tends to infinity: n — > oo. In such limit sum in 
formula (3.4) is replaced by integral over 3-dimensional space: 

(3.5) E{r)=Jp{r)^^d'r. 

Here p(f) is spatial density of charge at the point f . This value 
designates the amount of charge per unit volume. 

In order to find force acting on test charge q we should invert 
formula (3.3). As a result we obtain 

(3.6) F = gE(r). 

Force acting on a charge q in electric field is equal to the product 

of the quantity of this charge by the vector of intensity of field at 
the point, where charge is placed. However, charge q also produces 
electric field. Does it experience the action of its own field ? For 
point charges the answer to this question is negative. This fact 
should be treated as a supplement to principle of superposition. 
Total force acting on a system of distributed charges in electric 
field is determined by the following integral: 

(3.7) F = J p{r) E(r) d^r. 

Field E(r) in (3.7) is external field produced by external charges. 
Field of charges with density p{r) is not included into E(r). 

Concluding this section, note that formulas (3.4) and (3.5) hold 
only for charges at rest, which stayed at rest for sufficiently long 
time so that process of interaction transmitting reached the point 
of observation r. Fields produced by such systems of charges are 
called static fields, while branch of theory of electromagnetism 
studying such fields is called electrostatics. 
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j4. Lorentz force and Biot-Savart-Laplace law. 




Ampere law of interaction of parallel conductors with currents 
is an analog of Coulomb law for magnetic interactions. Accord- 
ing to near action principle, force F 
arises as a result of action of mag- 
netic field produced by a current in 
first conductor upon second conduc- 
tor. However, parallel conductors 
cannot be treated as point objects: 
formula (1.4) holds only for Z S> r. 
In order to get quantitative measure 
of magnetic field at some point r 
let's consider current I2 in (1.4) as 
a flow of charged particles of charge 
q each, and each moving along con- 
ductor with constant velocity v. If we denote by v the number 
of such particles per unit length of conductor, then in the whole 
length I we would have N = vl particles. Then during time in- 
terval t we would have n = vvt particles passing through a fixed 
cross-section of the conductor. They carry charge amounting to 
Q = qvvt. Therefore for current I2 in second conductor we get 

h = Q/t = qw. 

Upon calculating force acting on a segment of conductor of the 
length / by formula (1.4) we should divide it by the number of 
particles N contained in this segment. Then for the force per each 
particle we derive 



Fig. 4.1 



(4.1) 



2 hhl 2 hqv 



rN 



Formula determines (4.1) qualitative dependence of F on g and 
on v: each charged particle moving in magnetic field experiences 
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a force proportional to its charge q and to the magnitude of its 
velocity vector v = |v|, i.e. we have proportionality 



Force and velocity both are vectorial quantities. Simplest way to 
relate two vectors F and v is to use vector product of v with 
some third vectorial quantity H: 



Here c is scalar constant equal to light velocity. Vectorial quantity 
H(r) is a quantitative measure of magnetic field at the point r. 
It is called intensity of magnetic field at that point. Scalar factor 
1/c in (4.3) is used for to make H to be measured by the same 
units as intensity of electric field E in (3.6). Force F acting on 
a point charge in magnetic field is called Lorentz force. Total 
Lorentz force acting on a charge in electromagnetic field is a sum 
of two components: electric component and magnetic component: 



Formula (4.4) extends formula (3.6) for the case of general electro- 
magnetic fields. It holds not only for static but for time-dependent 
(non-static) fields. Surely the above derivation of formula (4.4) is 
empiric. Actually, one should treat formula (4.4) as experimental 
fact that do not contradict to another experimental fact (1.4) 
within theory being developed. 

Let's turn back to our conductors. Formula (4.3) can be 
interpreted in terms of currents. Each segment of unit length of 
a conductor with current I in magnetic field H experiences the 
force 



(4.2) 



F ~ qv. 



(4.3) 



F = ^[v, H(r)]. 



(4.4) 



F = gE+^[v, H]. 



(4.5) 
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acting on it. Here r is unit vector tangent to conductor and di- 
rected along current in it. Total force acting on circular conductor 
with current I is determined by contour integral 

(4.6) ^ = / ^ t"^^^^' "^""^^^^^ 

where s is natural parameter on contour (length) and r(s) is 
vector-function determining shape of contour in parametric form. 

Let's consider the case of two parallel conductors. Force F now 
can be calculated by formula (4.5) assuming that first conductor 
produces magnetic field H(r) that acts upon second conductor. 
Auxiliary experiment shows that vector H is perpendicular to 
the plane of these two parallel conductors. The magnitude of 
magnetic field if = |H| can be determined by formula (4.1): 

(4.7) H=^^-^. 

c r 

Here r is the distance from observation point to the conductor 
producing field at that point. 

Magnetic field produced by conductor with current satisfies 
superposition principle. In particular, field of infinite straight 
line conductor (4.7) is composed by fields produced by separate 
segments of this conductor. One cannot measure magnetic field of 
separate segment experimentally since one cannot keep constant 
current is such separate segment for sufficiently long time. But 
theoretically one can consider infinitesimally small segment of 
conductor with current of the length ds. And one can write 
formula for magnetic field produced by such segment of conductor: 

(4.8) rfH(r) = - ~ ''^ ds. 
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Here r is unit vector determining spatial orientation of infinitesi- 
mal conductor. It is always taken to be directed along current /. 
In practice, when calculating magnetic fields produced by circular 
conductors, formula (4.8) is taken in integral form: 



Like in (4.6), here s is natural parameter on the contour and f (s) 
is vectorial function determining shape of this contour. Therefore 
r(s) = dr{s)/ds. The relationship (4.8) and its integral form (4.9) 
constitute Biot-Savart-Laplace law for circular conductors with 
current. 

Biot-Savart-Laplace law in form (4.8) cannot be tested exper- 
imentally. However, in integral form (4.9) for each particular 
conductor it yields some particular expression for H(r). This 
expression then can be verified in experiment. 

Exercise 4.1. Using relationships (4.6) and (4.9), derive the 
law of interaction of parallel conductors with current in form (1.4). 

Exercise 4.2. Find magnetic field of the conductor with cur- 
rent having the shape of circle of the radius a. 



Conductors that we have considered above are kind of ideal- 
ization. They are linear, we assume them having no thickness. 
Real conductor always has some thickness. This fact is ignored 
when we consider long conductors like wire. However, in some 
cases thickness of a conductor cannot be ignored. For example, if 
we consider current in electrolytic bath or if we study current in 
plasma in upper layers of atmosphere. Current in bulk conductors 
can be distributed non-uniformly within volume of conductor. 
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§ 5. Current density 
and the law of charge conservation. 
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The concept of current density j is best one for describing such 
situation. 

Current density is vectorial quantity depending on a point of 
conducting medium: j = j(r). Vector of current density j(r) 
indicate the direction of charge transport at the point r. Its 
magnitude j = | j | is determined by the amount of charge passing 
through unit area perpendicular to vector j per unit time. Let's 
mark mentally some restricted domain within bulk conducting 
medium. Its boundary is smooth closed surface. Due to the 
above definition of current density the amount of charge flowing 
out from marked domain per unit time is determined by surface 
integral over the boundary of this domain, while charge enclosed 
within this domain is given by spatial integral: 



Here n is unit vector of external normal to the surface dCl 
restricting domain f2. 

Charge conservation law is one more fundamental experimental 
fact reflecting the nature of clcctromagnetism. In its classical 
form it states that charges cannot appear from nowhere and 
cannot disappear as well, they can only move from one point to 
another. Modern physics insert some correction to this statement: 
charges appear and can disappear in processes of creation and 
annihilation of pairs of elementary particles consisting of particle 
and corresponding antiparticlc. However, even in such creation- 
annihilation processes total balance of charge is preserved since 
total charge of a pair consisting of particle and antiparticle is 
always equal to zero. When applied to integrals (5.1) charge 
conservation law yields: Q = —J. This relationship means that 
decrease of charge enclosed within domain O is always due to 



(5.1) 
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charge leakage through the boundary and conversely increase of 
charge is due to incoming flow through the boundary of this 
domain. Let's write charge conservation law in the following 
form: 

(5.2) j^(^J pdh^+ J{i,n)dS = 0. 

n an 

Current density j is a vector depending on a point of conducting 
medium. Such objects in differential geometry arc called vector 
fields. Electric field E and magnetic field H are other examples 
of vector fields. Surface integral J in (5.1) is called flow of vector 
field j through the surface dfl. For smooth vector field any surface 
integral like J can be transformed to spatial integral by means of 
Ostrogradsky-Gauss formula. When applied to (5.2), this yields 

(5.3) J (§+divj) d'r = 0. 

n 

Note that Q in (5.3) is an arbitrary domain that we marked men- 
tally within conducting medium. This means that the expression 
being integrated in (5.3) should be identically zero: 

(5.4) ^ + div j = 0. 

The relationships (5.2) and (5.4) are integral and differential 
forms of charge conservation law respectively. The relationship 
(5.4) also is known as continuity equation for electric charge. 

When applied to bulk conductors with distributed current j 
within them, formula (4.6) is rewritten as follows: 



(5.5) 



F = /^[j(r), H(r)]d3r. 
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Biot-Savart-Laplace law for such conductors also is written in 
terms of spatial integral in the following form: 



In order to derive formulas (5.5) and (5.6) from formulas (4.6) 
and (4.8) one should represent bulk conductor as a union of linear 
conductors, then use superposition principle and pass to the limit 
by the number of linear conductors n — > oo. 



Let's consider some configuration of distributed charge with 
density p(r) which is concentrated within some restricted domain 
n. Let R be maximal linear size of the domain Let's choose 
coordinates with origin within this domain O and let's choose 
observation point r which is far enough from the domain of 
charge concentration: |r| ^ i?. In order to find electric field E(r) 
produced by charges in fl we use formula (3.5): 



Since domain Q in (6.1) is restricted, we have inequality |f| < R. 
Using this inequality along with |r| R, we can write Taylor 
expansion for the fraction in the expression under integration in 
(6.1). As a result we get power series in powers of ratio f/|r|: 



(5.6) 




§ 6. Electric dipole moment. 



(6.1) 




(6.2) 



r 



r 



r 



+ 



1 



WAlrl'M/ |r| 



r / r f \ f 



) 



+ 



|r — f |3 





Substituting (6.2) into (6.1), we get the following expression for 
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the vector of electric field E(r) produced by charges in O: 

. X , . r 3(r, D)r- IrPD 

(6.3) E(r) = Q^+ ^ j^|5 + • • • • 

First summand in (6.3) is Coulomb field of point charge Q placed 
at the origin, where Q is total charge enclosed in the domain O. 
It is given by integral (5.1). 

Second summand in (6.3) is known as field of point dipole 
placed at the origin. Vector D there is called dipole w.oment. For 
charges enclosed within domain O it is given by integral 

(6.4) ° " / '^(^)^'^^''- 

n 

For point charges dipole moment is determined by sum 

n 

(6.5) D = ^Qifi. 

i=l 

For the system of charges concentrated near origin, which is 
electrically neutral in whole, the field of point dipole 

, s . X 3(r, D)r- IrPD 

(6.6) E(r) = ^ ^^^3 ' ' 

is leading term in asymptotics for electrostatic field (3.4) or 
(3.5) as r — + oo. Note that for the system of charges with 
Q = dipole moment D calculated by formulas (6.4) and (6.5) 
is invariant quantity. This quantity remains unchanged when 
we move all charges to the same distance at the same direction 
without changing their mutual orientation: f — > f + tq. 
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Exercise 6.1. Concept of charge density is applicable to point 
charges as well. However, in this case p(r) is not ordinary function. 
It is distribution. For example point charge Q placed at the point 
r = is represented by density p{r) = Q 6{r), where S{r) is Dirac's 
delta-function. Consider the density 

(6.7) p(r) = <D, grad<5(r)) = 

i=l 

Applying formula (5.1), calculate total charge Q corresponding to 
this density (6.7). Using formula (6.4) calculate dipole moment 
for distributed charge (6.7) and Rnd electrostatic field produced 
by this charge. Compare the expression obtained with (6.6) and 
explain why system of charges described by the above density (6.7) 
is called point dipole. 

Exercise 6.2. Using formula (3.7) find the force acting on 
point dipole in external electric field E(r). 

§ 7. Magnetic moment. 

Let's consider situation similar to that of previous section. 
Suppose some distributed system of currents is concentrated in 
some restricted domain near origin. Let R be maximal linear 
size of this domain f2. Current density j(r) is smooth vector- 
function, which is nonzero only within f2 and which vanishes at 
the boundary dQ and in outer space. Current density j(r) is 
assumed to be stationary, i. e. it doesn't depend on time, and it 
doesn't break charge balance, i. e. p{r) = 0. Charge conservation 
law applied to this situation yields 



(7.1) 



divj = 0. 
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In order to calculate magnetic field H(r) we use Biot-Savart- 
Laplace law written in integral form (5.6): 

(7.2) H(r)= /iiiilll^^d^f. 



n 



Assuming that |r| » R, we take Taylor expansion (6.2) and 
substitute it into (7.2). As a result we get 



(7 3) 

^ /■ 3(r,f)[j(r),r]-|rnj(r),f] ^3. ^ 

J c|r|5 
n 

Lemma 7.1. First integral in (7.3) is identically equal to zero. 

Proof. Denote this integral by Hi(r). Let's choose some 
arbitrary constant vector e and consider scalar product 

(7.4) (H„ e) = / ^ / (m, [r e]) ^3,_ 

Q n 

Then define vector a and function /(f) as follows: 

[r, e] / 
a=TM3' /(r) = (a, r). 



Vector a does not depend on f , therefore in calculating integral 
(7.4) we can take it for constant vector. For this vector we derive 
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a = grad/. Substituting this formula into the (7.4), we get 

<Hi,e) = J{j, gr&d f)d'i = 
n 

= J div(/ j) _ y /divjd^f. 



(7.5) " 



Last integral in (7.5) is equal to zero due to (7.1). Previous inte- 
gral is transformed to surface integral by means of Ostrogradsky- 
Gauss formula. It is also equal to zero since j(f) vanishes at the 
boundary of domain Cl. Therefore 

(7.6) (Hi,e) = j f{j,n)dS = 0. 

on 

Now vanishing of vector Hi(r) follows from formula (7.6) since e 
is arbitrary constant vector. Lemma 7.1 is proved. □ 

Let's transform second integral in (7.3). First of all we denote 
it by H2(r). Then, taking an arbitrary constant vector e, we form 
scalar product (H2, e). This scalar product can be brought to 

(7.7) (H2,e> = ^y(j(f),b(f))d3f, 



where b(f) = 3(r, f) [r, e] — |rp [f, e]. If one adds gradient of 
arbitrary function /(f) to b(f), this wouldn't change integral 
in (7.7). Formulas (7.5) and (7.6) form an example of such 
invariance. Let's specify function /(f), choosing it as follows: 

(7.8) /(f) = -^(r,f)(f, [r,e]). 
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For gradient of function (7.8) by direct calculations we find 

grad/(f) = -^(f, [r,e])r- ^(r, f> [r,e] = 

= -3 (r, f ) [r, e] - ^ (r (f , [r, e]) - [r, e] (r, f )) . 

Now let's use well-known identity [a, [b, c]] = b (a, c) — c (a, b). 
Assuming that a = f , b = r, and c = [r, e], we transform the 
above expression for grad / to the following form: 

(7.9) grad /(f) = -3 (r, f) [r,e] - ^ [f, [r, [r, e]]]. 

Right hand side of (7.9) contains triple vectorial product. In order 
to transform it we use the identity [a, [b, c]] = b (a, c) — c (a, b) 
again, now assuming that a = r, b = r, and c = e: 

grad /(f) = -3 (r, f) [r,e] - ^ (r,e) [f, r] + ^ \rf [f, e]. 

Let's add this expression for grad/ to vector b(f). Here is 
resulting new expression for this vector: 

(7.10) b(f) = -^(r,e)[f,r] + i|rnf,e]. 
Let's substitute (7.10) into formula (7.7). This yields 

-3(r,e)(r, [j(f), f]) + |rp (e, [ J (?),?]) 3- 

2c|r|5 
a 

Note that quantities j(f ) and f enter into this formula in form of 



(H2, e) = / 
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vector product [j(f), f] only. Denote by M the following integral: 
(7.11) M = yi^^i^d3~_ 



2c 



Vector M given by integral (7.11) is called magnetic moment for 
currents with density j(f). In terms of M the above relationship 
for scalar product (H2, e) is written as follows: 

(7^12) (H..e> = ='<'-°><'-">-l'l' 



If we remember that e in formula (7.12) is an arbitrary constant 
vector, then from (7.3) and lemma 7.1 we can conclude that the 
field of point magnetic dipole 

, . 3(r, M)r- IrPM 
(7-13) H(r) = ^ ' ' ' 

is leading term in asymptotical expansion of static magnetic field 
(4.9) and (5.6) as r — *• 00. 

Like electric dipole moment D of the system with zero total 
charge (5 = 0, magnetic moment M is invariant with respect 
to displacements r ^ r + Tq that don't change configuration of 
currents within system. Indeed, under such displacement integral 
(7.11) is incremented by 



(7.14) AM = j 



2c 



(fr = 0. 



Integral in formula (7.14) is equal to zero by the same reasons as 
in proof of lemma 7. 1 . 
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Exercise 7.1. Consider localized system of currents ]{r) with 
current density given by the following distribution: 



Verify the relationship (7.1) for the system of currents (7.15) and 
find its magnetic moment M. Applying formula (5.6), calculate 
magnetic field of this system of currents and explain why this sys- 
tem of currents is called point magnetic dipole. 

Exercise 7.2. Using formula (5.5), find the force acting upon 
point magnetic dipole in external magnetic field H(r). 

Exercise 7.3. By means of the following formula for the torque 



find torque M. acting upon point magnetic dipole (7.15) in homo- 
geneous magnetic held H = const. 



Remember that we introduced the concept of flow of vector 
field through a surface in considering charge conservation law (see 
integral J in (5.1)). Now we consider flows of vector flelds E(r) 
and H(r), i.e. for electric field and magnetic field: 



(7.15) 



j(r) = -c[M, grad5(r)]. 




§ 8. Integral equations 
for static electromagnetic field. 



(8.1) 





s 



s 



Let S be closed surface enveloping some domain fl, i.e. S = 
dQ. Electrostatic field E is determined by formula (3.5). Let's 
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substitute (3.5) into first integral (8.1) and then let's change order 
of integration in resulting double integral: 



12) 



dn 



Inner surface integral in (8.2) is an integral of explicit function. 
This integral can be calculated explicitly: 



.3) / 



(r-f, nW) ^^_ j 0, iff ^O, 
|r - f |3 I 47r, if f e n. 

dQ ^ 

Here by f2 = i7 U dfl we denote closure of the domain O. 

In order to prove the relationship (8.3) let's consider vector 
field m(r) given by the following formula: 

(8.4) m(r) = 



Vector field m(r) is smooth everywhere except for one special 
point r = f . In all regular points of this vector field by direct 
calculations we find divm = 0. If f ^ 57 special point of the field 
m is out of the domain ^l. Therefore in this case we can apply 
Ostrogradsky- Gauss formula to (8.3): 

j{ni, n)dS = y^divmrf^r = 0. 

dn n 

This proves first case in formula (8.3). In order to prove second 

case, when f G 51, we use tactical maneuver. Let's consider spher- 
ical e-neighborhood O = of special point r = f . For sufficiently 
small e this neighborhood O = is completely enclosed into the 
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domain Q. Then from zero divergency condition divm = for 
the field given by formula (8.4) we derive 

(8.5) j{m,n)dS=j{ni,n)dS = iTr. 

on do 

The value of last integral over sphere dO in (8.5) is found by direct 
calculation, which is not difficult. Thus, formula (8.3) is proved. 
Substituting (8.3) into (8.2) we get the following relationship: 



j (E, n)dS = 4T: j p{r)d^v. 



(8.6) 

dn 

This relationship (8.6) can be formulated as a theorem. 

Theorem (on the flow of electric field). Flow of electric Held 
through the boundary of restricted domain is equal to total charge 
enclosed within this domain multiplied by A-k. 

Now let's consider flow of magnetic field Ti in (8.1). Static 
magnetic field is determined by formula (5.6). Let's substitute 
H(r) given by (5.6) into second integral (8.1), then change the 
order of integration in resulting double integral: 

,,,, „.//i W.;:;};-w> .s.3, 

on 

It's clear that in calculating inner integral over the surface 
vector j can be taken for constant. Now consider the field 

(8.8) m(r) = 

c r — r 
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Like (8.4), this vector field (8.8) has only one singular point r = f . 
Divergency of this field is equal to zero, this fact can be verified 
by direct calculations. As appears in this case, singular point 
makes no effect to the value of surface integral in (8.7). Instead 
of (8.3) in this case we have the following formula: 

(8.9) n ([j^^-^W'-)) ,s=0. 

J c |r-r|'^ 
an 

For r ^ Q the relationship (8.9) follows from divm = by 
applying Ostrogradsky-Gauss formula. For f G f7 we have the 
relationship similar to the above relationship (8.5): 

(8.10) j (m, n) dS = j (m, n) dS = 0. 

an do 

However, the value of surface integral over sphere dO in this case 
is equal to zero since vector m(r) is perpendicular to normal 
vector n at all points of sphere dO. As a result of substituting 
(8.9) into (8.7) we get the relationship 

(8.11) j{U,n)dS = 0, 

an 

which is formulated as the following theorem. 

Theorem (on the flow of magnetic field). Total Bow of mag- 
netic field through the boundary of any restricted domain is equal 
to zero. 

Let r(s) be vectorial parametric equation of some closed spatial 
curve r being the rim for some open surface S, i. e. F = dS. Open 
surface S means that S and F have empty intersection. By S we 
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denote the closure of the surface S. Then S = S U T. Taking 
s for natural parameter on F, we define circulation for electric 

and magnetic fields in form of the 
following contour integrals: 




i.U) 



e = ^(E, t) ds, 
r 

() = ^(H, r) ds. 



Substituting (3.5) into (8.12) and 
Fig. 8.1 changing the order of integration in 

resulting double integral, we get the 
following equality for circulation of electric field: 



(8.13) 



p(f) 



(r(s) -f, t(s)) 



\v{s) 



f|3 



dsd^Y. 



Due to (8.13) wc need to consider vector field (8.4) again. For 
f ^ F, taking into account T = dS and applying Stokes formula, 
we can transform contour integral in (8.13) to surface integral: 



(8.14) 



(r(5) -r, t{s)) 
|r(s) — f |3 



ds 



j (rotm, n)^^ = 0. 



Values of integral (8.14) at those points f £ F arc of no matter 
since when substituting (8.14) into integral (8.13) such points 
constitute a set of zero measure. 

Vanishing of integral (8.14) for f ^ F follows from rotm = 0, 
this equality can be verified by direct calculations. Singular point 
r = f of vector field (8.4) is unessential since surface S, for which 
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r is a boundary, can be deformed so that i ^ S. The result of 
substituting (8.14) into (8.13) can be written as an equation: 



(8.15) 




dS 



Theorem (on the circulation of electric field). TotaJ circula- 
tion of static electric field along the boundary of any restricted 
open surface is equal to zero. 

Formula like (8.15) is available for magnetic field as well. Here 
is this formula that determines circulation of magnetic field: 

(8.16) ^(H, T)ds=^ J ( j, n) dS. 

ds s 

Corresponding theorem is stated as follows. 

Theorem (on the circulation of magnetic field). Circulation of 

static magnetic Geld along boundary of restricted open surface is 
equal to total electric current penetrating this surface multiplied 
by fraction 4 tt/c. 

Integral over the surface 5 now is in right hand side of formula 

(8.16) explicitly. Therefore surface spanned over the contour T 
now is fixed. We cannot deform this surface as we did above 
in proving theorem on circulation of electric field. This leads to 
some technical complication of the proof. Let's consider e-blow- 
up of surface S. This is domain ^{s) being union of all £-balls 
surrounding all point r e 5. This domain encloses surface S and 
contour F = dS. If e ^ 0, domain il{e) contracts to S. 

Denote by D{£) = M.^ \ ^{e) exterior of the domain ^{e) and 
then consider the following modification of formula (5.6) that 
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expresses Biot-Savart-Laplace law for magnetic field: 

(8.17) H(r) = lim /' ^Ml^^^d^f. 
^ ' ^ ' E^o J c |r - r|3 

D{e) 

Let's substitute (8.17) into integral (8.12) and change the order of 
integration in resulting double integral. As a result we get 

(8.18) f, = lin. / /l([j(^)^;(^)-^^];^(^)),,,3,. 
' e^o J J c \r{s) - r|3 

Die) r 

In inner integral in (8.18) we see vector field (8.8). Unlike vector 
filed (8.4), rotor of this field is nonzero: 

3 (r — f , j) (r — f) — |r — f P j 

(8.19) rotm = — ^ ' ^ ^ ' ' 



Using Stokes formula and taking into account (8.19), we can 
transform contour integral (8.18) to surface integral: 

l ([j(f), r(.s)-r],T(5)) ^_ 
c |r(s) — f |3 

3 (r - f, j(f)) (r - f , n(r)) - |r - f ^ ( j(f), n(r)) 



/ 



■dS. 

c |r — r|° 

s 



Denote by m(f ) vector field of the following form: 

„ 3(f -r, n(r)) (f -r) - |f -r|2n(r) 
m(r) = — ^ 
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In terms of the field m(f ) formula for \) is written as 
i) = l^^j y(m(f)J(f)>d5d3f. 

D(e) S 

Vector field rn(f) in tliis formnla lias cubic singularity |f — r|~"^ 
at the point f = r. Such singularity is not integrable in M.^ (if 
we integrate with respect to d^r). This is why we use auxiliary 
domain D{e) and limit as £ ^ 0. 

Let's change the order of integration in resulting double integral 
for circulation f). This leads to formula 

J J <m(f)J(f))d3fd5 = J J (grad/(f),j(f))d3fd5, 

S D{e) S D{s) 

since vector field m(f) apparently is gradient of the function /(f): 

(f - r, n(r)> 



(8.20) /(f) 



c r ■ 



Function /(f) vanishes as f — > oo. Assume that current density 
also vanishes as f — > oo. Then due to the same considerations as 
in proof of lemma 7.1 and due to formula (7.1) spatial integral in 
the above formula can be transformed to surface integral: 

(8.21) '5 = 1™/ / /(f)(j(f), ii(f))rf5d5. 

S dD(e) 



Let's change the order of integration in (8.21) then take into 
account common boundary dD{e) = dCl{e). Outer normal to the 
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surface dD{e) coincides with inner normal to 9r2(e). This coin- 
cidence and expUcit form of ftmction (8.20) lead to the following 
expression for circulation of magnetic field f) : 

f (j(f), n(f)) /"(f-r, n(r)) 
(8.22 f) = lim / ^-"^ ^' ^ / ^ ' dSdS. 
' e^o J c y |r-r|3 

an(e) s 
Let's denote by V{t) inner integral in formula (8.22): 

s 

Integral (8.23) is well-known in mathematical physics. It is called 
potential of double layer. There is the following lemma, proof of 
which can be found in [1]. 

Lemma 8.1. Double layer potential (8.23) is restricted func- 
tion in R3 \^ s_ yi^t eacJi inner point r G S there are side limits 

14(f) = lim V(t), 

inner limit V-{r) as r tends to r G S from inside along normal 
vector n, and outer limit V+ (f ) as r tends to r € S from outside 
against the direction of normal vector n. Thereby V+ — V- = 
for all points r G S. 

In order to calculate limit in formula (8.22) we need to study 
the geometry of £-blow-up of the surface S. On Fig. 8.2 below 
we see cross-section of the domain Cl{e) obtained from the surface 
S shown on Fig. 8.1. For sufficiently small s boundary of the 
domain Q{e) is composed of three parts: 



(8.24) 



dfl{s) = SoUS+US- 
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Surface Sq is a part of e-blow-up of the contour T Area of this 
surface Sq satisfies the relationship 



(8.25) 




So ~ swL as s 
S+ 



0, 



Fig. 8.2 



where L is length of contour 
r. Surfaces 5+ and S- are ob- 
tained as a result of normal shift 
of surface 5* to the distance e 
along normal vector n, and to 
the same distance against nor- 
mal vector n. 

Substituting (8.24) into (8.22) 
we break surface integral over 
9f2(e) into three parts. Since 
double layer potential and func- 
tion |j(f)| are restricted, we get 
the relationship 



(8.26) 



lim 



V{i) 



dS = Q. 



So 



For other two summand we also can calculate limits as e — > 0: 



(8.27) 



S± 



dS 



: J 14 (r) 



(j(r), n(r)) 



dS. 



We shall not load reader with the proof of formulas (8.24), (8.25) 
and (8.27), which are sufficiently obvious. Summarizing (8.26) 
and (8.27) and taking into account lemma 8.1, we obtain 



(8.28) 



f) = f /<j(r),n(r))d5. 
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This relationship (8.28) completes derivation of formula (8.16) 
and proof of theorem on circulation of magnetic field in whole. 

Exercise 8.1. Verify the relationship divm = for vector 
fields (8.4) and (8.8). 

Exercise 8.2. Verify the relationship (8.19) for vector field 
given by formula (8.8). 

Exercise 8.3. Calculate grad/ for the function (8.20). 



In §8 we have derived four integral equations for electric and 
magnetic fields. They are used to be grouped into two pairs. 
Equations in first pair have zero right hand sides: 



Right hand sides of equations in second pair are non-zero. They 
are determined by charges and currents: 



§ 9. Differential equations 
for static electromagnetic field. 



(9.1) 




dS 





an 



n 



(9.2) 




OS 



s 



Applying Ostrogradsky-Gauss formula and Stokes formula, one 
can transform surface integrals to spatial ones, and contour inte- 
grals to surface integrals. Then, since O is arbitrary domain and 
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S is arbitrary open surface, integral equations (9.1) and (9.2) can 
be transformed to differential equations: 

(9.3) divH = 0, rotE = 0, 

47r 

(9.4) div E = 4ttp, rot H = — j. 

When considering differential equations (9.3) and (9.4), we should 
add conditions for charges and currents being stationary: 

(-) 

The relationship (7.1) then is a consequence of (9.5) and charge 

conservation law. 

Differential equations (9.3) and (9.4) form complete system 
of differential equations for describing stationary electromagnetic 
fields. When solving them functions p(r) and j(r) arc assumed to 
be known. If they are not known, one should have some additional 
equations relating p and j with E and H. These additional 
equations describe properties of medium (for instance, continuous 
conducting medium is described by the equation j = a E, where a 
is conductivity of medium). 
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CHAPTER II 



CLASSICAL ELECTRODYNAMICS 



§ 1. Maxwell equations. 



Differential equations (9.3) and (9.4), which we have derived 
in the end of Chapter I, describe fields generated by stationary 
charges and currents. They are absolutely unsuitable if we are 
going to describe the process of haw electromagnetic interaction is 
transmitted in space. Note that the notion of field was introduced 
within framework of the concept of near action for describing 
the object that transmit interaction of charges and currents. For 
static fields this property is revealed in a very restrictive form, 
i. e. we use fields only to divide interaction of charges and currents 
into two processes: creation of a field by charges and currents is 
first process, action of this field upon other currents and charges 
is second process. Dynamic properties of the field itself appears 
beyond our consideration. 

More exact equations describing process of transmitting elec- 
tromagnetic interaction in its time evolution were suggested by 
Maxwell. They are the following ones: 



(1.2) 



(1.1) 



divH = 0, 



div E = 47rp, 



rotH 



rotE 




1 dU 
4tt . 1 dE 
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It is easy to see that equations (1.1) and (1.2) are generalizations 
for the (9.3) and (9.4) from Chapter I. They are obtained from 
latter ones by modifying right hand sides. Like equations (9.3) 
and (9.4) in Chapter I, Maxwell equations (1.1) and (1.2) can be 
written in form of integral equations: 



(1.3) 



(1.4) 



j (H, n) dS = 0, 
an 

/(E, T)d.= -i||(H, n)dS, 
as s 

j (E, n) dS = 4TT j pd^r, 
an n 

/(H, T)ds = ^l(j,n)dS+\fJ(^,n) dS. 



Consider contour integral in second equation (1.3). Similar 
contour integral is present in second equation (1.4). However, 
unlike circulation of magnetic field, circulation of electric field 




as 



possess its own physical interpretation. If imaginary contour 
r = dS in space is replaced by real circular conductor, then 
electric field with nonzero circulation induces electric current in 
conductor. The quantity e from (1.5) in this case is called electro- 
motive force of the field E in contour. Electromotive force e ^ 
in contour produce the same effect as linking electric cell with 
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voltage e into this contour. Experimentally it reveals as follow: 
alternating magnetic field produces electric field with nonzero cir- 
culation, this induces electric current in circular conductor. This 
phenomenon is known as electromagnetic induction. It was first 
discovered by Faraday. Faraday gave qTialitative description of 
this phenomenon in form of the following induction law. 

Faraday's law of electromagnetic induction. Electromo- 
tive force of induction in circular conductor is proportional to the 
rate of changing of magnetic flow embraced by this conductor. 

Faraday's induction law was a hint for Maxwell when choosing 
right hand side in second equation (1.1). As for similar term in 
right hand side of second equation (1.2), Maxwell had written it 
by analogy. Experiments and further development of technology 
proved correctness of Maxwell equations. 

Note that charge conservation law in form of relationship (5.4) 
from Chapter I is a consequence of Maxwell equations. One 
should calculate divergency of both sides of second equation (1.2): 

„ 47r ,. . 1 adivE 

divrotH = — divj H — — , 

c c at 

then one should apply the identity divrotH = 0. When combined 
with the first equation (1.2) this yields exactly the relationship 
(5.4) from Chapter I. 

Equations (1.1) and (1.2) form complete system for describ- 
ing arbitrary electromagnetic fields. In solving them functions 
p{r,t) and j{r,t) should be given, or they should be determined 
from medium equations. Then each problem of electrodynam- 
ics mathematically reduces to some boundary-value problem or 
mixed initial- value/boundary- value problem for Maxwell equa- 
tions optionally completed by medium equations. In this section 
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we consider only some very special ones among such problems. 
Our main goal is to derive important mathematical consequences 
from Maxwell equations and to interpret their physical nature. 



Suppose that in bulk conductor we have a current with density 
j, and suppose that this current is produced by the flow of charged 
particles with charge q. If z/ is the number of such particles per 
unit volume and if v is their velocity, then j = qw. Recall that 
current density is a charge passing through unit area per unit 
time (see § 5 in Chapter I) . 

In electromagnetic field each particle experiences Lorentz force 
determined by formula (4.4) from Chapter I. Work of this force 
per unit time is equal to (F, v) = g (E, v). Total work produced 
by electromagnetic field per unit volume is obtained if one multi- 
plies this quantity by i^, then w = qv (E, v) = (E, j). This work 
increases kinetic energy of particles (particles are accelerated by 
field). Otherwise this work is used for to compensate forces of 
viscous friction that resist motion of particles. In either case 
total power spent by electromagnetic field within domain O is 
determined by the following integral: 



Let's transform integral (2.1). Let's express current density j 
through E and H using second equation (1.2) for this purpose: 



§ 2. Density of energy and energy flow 
for electromagnetic field. 



(2.1) 




n 



(2.2) 
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Substituting this expression (2.2) into formula (2.1), we get 

(2.3) = ^/(E, rotH)d3r - -L J |(e, E) d^r. 

n n 

In order to implement further transformations in formula (2.3) 
we use well-known identity div [a, b] = (b, rota) — (a, rotb). 
Assuming a = H and b = E, for we get 

W=^ [ div[H, E] d^r+^ [ (H, rot E> d^r - A / ^ d^r. 
Att J ^ ' ^ in J ^ ' ' dtJ 8tt 

Q n Q 

First integral in this expression can be transformed by means of 
Ostrogradsky-Gauss formula, while for transforming rotE one can 
use Maxwell equations (1.1): 

M VV- + /^([E.H|,„>.S + |/E±m#r = 0. 

an n 

Let's denote by S and e vectorial field and scalar field of the form 
(2.5, S=flE.H,, 

The quantity e in (2.5) is called density of energy of electromag- 
netic field. Vector S is known as density of energy flow. It also 
called Umov- Pointing vector. Under such interpretation of quan- 
tities (2.5) the relationship (2.4) can be treated as the equation 
of energy balance. First summand in (2.4) is called dissipation 
power, this is the amount of energy dissipated per unit time at 
the expense of transmitting it to moving charges. Second sum- 
mand is the amount of energy that flows from within domain O to 
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outer space per unit time. These two forms of energy losses lead 
to diminishing the energy stored by electromagnetic field itself 
within domain (see third summand in (2.4)). 

Energy balance equation (2.4) can be rewritten in differential 
form, analogous to formula (5.4) from Chapter I: 

de 

(2.6) — +divS + «; = 0. 

at 

Here w =(E, j) is a density of energy dissipation. Note that in 
some cases w and integral (2.1) in whole can be negative. In such 
a case we have energy pumping into electromagnetic field. This 
energy then flows to outer space through boundary of the domain 
VL. This is the process of radiation of electromagnetic waves from 
the domain Q.. It is realized in antennas (aerials) of radio and 
TV transmitters. If we eliminate or restrict substantially the 
energy leakage from the domain f2 to outer space, then we would 
have the device like microwave oven, where electromagnetic field 
is used for transmitting energy from radiator to beefsteak. 

Electromagnetic field can store and transmit not only the 
energy, but the momentum as well. In order to derive momentum 
balance equations let's consider again the current with density j 
due to the particles with charge q which move with velocity v. 
Let V be concentration of these particles, i. e. number of particles 
per unit volume. Then '} = qw and p = qv. Total force acting 
on all particles within domain Q. is given by integral 

(2.7) F = J pEdh + J ^[i,n] 

n n 

In order to derive formula (2.7) one should multiply Lorentz force 
acting on each separate particle by the number of particles per 
unit volume u and then integrate over the domain Cl. 
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Force F determines the amount of momentum transmitted 
from electromagnetic field to particles enclosed within domain Q. 
Integral (2.7) is vectorial quantity. For further transformations of 
this integral let's choose some constant unit vector e and consider 
scalar product of this vector e and vector F: 

(2.8) <F, e) = Jp{ E, e) d'r + J ^ (e, [j, H]) d^r. 

n 

Substituting (2.2) into (2.8), we get 

(F, J p (E, e) d^r + ^ J (e, [rot H, H]) d^r - 
n n 

(2.9) 



Recalling well-known property of mixed product, we do cyclic 
transposition of multiplicands in second integral (2.9). Moreover, 
we use obvious identity [dE/dt, U] = d [E, H] /dt - [E, dU/dt]. 
This yields the following expression for (F, e): 



(F, e) = J p{ E, e) d\ + ^ j (rot H, [H, e]) rf^r - 

n n 
Now we apply second equation of the system (1.1) written as 
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dn/dt = — c rotE. Then we get formula 



(2.10) 



+ 



/ 



(rot H, [H, e]> + (rot E, [E, e]> 



47r 



In order to transform last two integrals in (2.10) we use the 
following three identities, two of which we already used earlier: 

[a, [b, c]] = b(a, c) - c (a, b), 
(2.11) div [a, b] = (b, rot a) - (a, rot b), 

rot [a, b] = a div b — b diva — {a, b}. 



Here by curly brackets wc denote commutator of two vector 
fields a and b (see [2]). Traditionally square brackets are used 
for commutator, but here by square brackets we denote vector 
product of two vectors. From second identity (2.11) we derive 

(rotH, [H, e]) = div [H, [H, e]] + (H, rot [H, e]). 

In order to transform rot[H, e] we use third identity (2.11): 
rot[H, e] = -e divH - {H, e}. Then 

(H, rot [H, e]) = -(H, e) divH + ^ ^ — = 
= -(H, e) divH+^(e, grad|H|2). 
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Let's combine two above relationships and apply first identity 
(2.11) for to transform double vectorial product [H, [H, e]] in 
first of them. As a result we obtain 

(rotH, [H, e]) =div(H(H, e)) - div(e |H|2) - 

-(H, e> divH + i(e, grad|H|2>. 

But div(e |Hp)= (e, grad|Hp). Hence as a final result we get 
(rotH, [H, e]) = -(H, e) divH + 

(2.12) 



+ div (u (H, e> - ^ e |H|2^ . 



Quite similar identity can be derived for electric field E: 
(rotE, [E, e]) = -(E, e) divE + 

(2.13) 



+ div(^E(E, e)-ie|E|2^. 



The only difference is that due to Maxwell equations divH = 0, 
while divergency of electric field E is nonzero: divE = Airp. 

Now, if we take into account (2.12) and (2.13), formula (2.10) 
can be transformed to the following one: 



an 



(E, e)(n, E) + (H, e)(n, H) ^ 



J Stt dt J A-KC 

Denote by a linear operator such that the result of applying this 
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operator to some arbitrary vector e is given by formula 

(2.14) ^^_E(E,e)+H(H,e) lEP + IHP 

^ ^ 47r Stt 

Formula (2.14) defines tensorial field a of type (1, 1) with the 
following components: 

(2.15) 

Tensor cr with components (2.15) is called tensor of the density 
of momentum flow. It is also known as Maxwell tensor. Now let's 
define vector of momentum density p by formula 

(2.16) p 



47rc 



In terms of the notations (2.15) and (2.16) the above relationship 
for (F, e) is rewritten as follows: 

(2.17) (F, e) + j{ae,u)dS+jJ (p, e> d^r = 0. 

Operator of the density of momentum flow cr is symmetric, i. e. 
(ere, n) = (e, an). Due to this property and because e is 
arbitrary vector we can rewrite (2.17) in vectorial form: 

(2.18) F + j andS + j pd^r = 0. 

This equation (2.18) is the equation of momentum balance for 
electromagnetic field. Force F, given by formula (2.7) determines 
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loss of momentum stored in electromagnetic field due to transmit- 
ting it to moving particles. Second term in (2.18) determines loss 
of momentum due to its flow through the boundary of the domain 
fl. These two losses lead to diminishing the momentum stored 
by electromagnetic field within domain Q (see third summand in 
(2.18)). 

The relationship (2.18) can be rewritten in differential form. 
For this purpose we should define vectorial divergency for tensorial 
field (T of the type (1, 1). Let 



Eda. 

i=l 



""j 



Then differential form of (2.18) is written as 

(2.20) ^+divcr + f = 0, 

where f = + [j, H]/c is a density of Lorentz force, while 
vectorial divergency is determined according to (2.19). 

Thus, electromagnetic field is capable to accumulate within 
itself the energy and momentum: 

(2.21) 8=[\-^l±&d^r, P=/M.3r. 
^ ' J J Attc 

n n 

It is also capable to transmit energy and momentum to material 
bodies. This confirms once more our assertion that electromag- 
netic field itself is a material entity. It is not pure mathematical 
abstraction convenient for describing interaction of charges and 

currents, but real physical object. 

Exercise 2.1. Verify that relationships (2.11) hold. Check on 
the derivation of {2.12) and (2.13). 
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§ 3. Vectorial and scalar potentials 
of electromagnetic field. 

In section 2 we have found that electromagnetic field possess 
energy and momentum (2.21). This is very important conse- 
quence of Maxwell equations (1.1) and (1.2). However we have 
not studied Maxwell equations themselves. This is system of four 
equations, two of them are scalar equations, other two are vecto- 
rial equations. So they are equivalent to eight scalar equations. 
However we have only six undetermined functions in them: three 
components of vector E and three components of vector H. So 
observe somewhat like excessivencss in Maxwell equations. 

One of the most popular ways for solving systems of algebraic 
equations is to express some variable through other ones by 
solving one of the equations in a system (usually most simple 
equation) and tlic;n substituting the expression obtained into other 
equations. Thus we exclude one variable and diminish the number 
of equations in a system also by one. Sometimes this trick 
is applicable to differential equations as well. Let's consider 
Maxwell equation div H = 0. Vector field with zero divergency is 
called vortex field. For vortex fields the following theorem holds 
(see proof in book [3]). 

Theorem on vortex field. Each vortex field is a rotor of some 
otlier vector field. 

Let's write the statement of this theorem as applied to mag- 
netic field. It is given by the following relationship; 

(3.1) H = rotA. 

Vector field A, whose existence is granted by the above theorem, 
is called vector-potential of electromagnetic field. 
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Let's substitute vector H as given by (3.1) into second Maxwell 
equation (1.1). This yields the equality 

(3.2) rotE+i-^rotA = rot( E+-^ ) =0. 

cot \ c at J 

Vector field with zero rotor is called potential fi,eld. It is vector 
field E + {dA/dt)/c in formula (3.2) which is obviously potential 
field. Potential fields are described by the following theorem (see 
proof in book [3]). 

Theorem on potential field. Each potential field is a gradi- 
ent of some scalar field. 

Applying this theorem to vector field (3.2), we get the relation- 
ship determining scalar potential of electromagnetic field (p: 

(3.3) E+-— = -grad<p. 

Combining (3.1) and (3.3), we can express electric and magnetic 
fields E and H through newly introduced fields A and ip: 

, l^A 

(3.4) ^=-S™^^-ca^' 
H = rot A. 

Upon substituting (3.4) into first pair of Maxwell equations 
(1.1) we find them to be identically fulfilled. As for second pair of 
Maxwell equations, substituting (3.4) into these equations, we get 

— Aip — div A = 4 TT p, 

^^■^^ ... * A* 15 , 1 S^A 47rj 

graddiv A - AA + -- grad^ + -_ = _. 
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In deriving (3.5) we used relationships 



(3.6) 



div grad <^ = A<^, 

rot rot A = grad div A — AA. 



Second order differential operator A is called Laplace operator. In 
rectangular Cartesian coordinates it is defined by formula 



In order to simplify the equations (3.5) we rearrange terms in 
them. As a result we get 



Differential equations (3.8) are Maxwell equations written in 
terms of A and ip. This is system of two equations one of 
which is scalar equation, while another is vectorial equation. As 
we can see, number of equations now is equal to the number of 
undetermined functions in them. 

§ 4. Gauge transformations and Lorentzian gauge. 

Vectorial and scalar potentials A and ip were introduced in 
§3 as a replacement for electric and magnetic fields E and H. 
However, fields A and (p arc not physical fields. Physical fields 
E and H are expressed through A and ip according to formulas 
(3.4), but backward correspondence is not unique, i.e. fields A 
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and if are not uniquely determined by physical fields E and H. 
Indeed, let's consider transformation 

A = A + gradi/', 
(4-1) 1 d^j 

c at 

where ip{''^jt) is an arbitrary function. Substituting (4.1) into 
formula (3.4), we immediately get 

E = E, H = H. 

This means that physical fields E, H determined by fields A, (f 
and by fields A, ip do coincide. Transformation (4.1) that do not 
change physical fields E and H is called gauge transformation. 

We use gauge transformations (4.1) for further simplification of 
Maxwell equations (3.8). Let's consider the quantity enclosed in 
brackets in right hand sides of the equations (3.8): 

Denote by □ the following differential operator; 

1 

Operator (4.3) is called d'Alambert operator or wave operator. 
Differential equation Du = v is called d'Alambert equation. 

Using gauge freedom provided by gauge transformation (4.1), 
we can fulfill the following condition: 

(4.4) -$+divA = 0. 

c ot 
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For this purpose we should choose V' solving d'Alambert equation 

°^ = -Gf +'^^^^)- 

It is known that d'Alambert equation is solvable under rather 
weak restrictions for its right hand side (see book [1]). Hence 
practically always we can fulfill the condition (4.4). This condition 
is called Lorentzian gauge. 

If Lorentzian gauge condition (4.4) is fulfilled, then Maxwell 
equations (3.8) simplify substantially: 



(4.5) □(p = 47rp, □A=^. 



They look like pair of independent d'Alambert equations. How- 
ever, one shouldn't think that variables A and <p are completely 
separated. Lorentzian gauge condition (4.4) itself is an additional 
equation requiring concordant choice of solutions for d'Alambert 
equations (4.5). 

D'Alambert operator (4.3) is a scalar operator, in (4.5) it acts 
upon each component of vector A separately. Therefore operator 
□ commutates with rotor operator and with time derivative as 
well. Therefore on the base of (3.4) we derive 



47r di 47r 

(4.6) DE = -4tt grad p ^ DH = — rot j. 

at c 



These equations (4.6) have no entries of potentials A and ip. 
They arc written in terms of real physical fields E and H, and 
are consequences of Maxwell equations (1.1) and (1.2). However, 
backward Maxwell equations do not follow from (4.6). 
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§ 5. Electromagnetic waves. 

In previous Chapter we considered static electromagnetic fields. 
Such fields are uniquely determined by static configuration of 

charges and currents (see formu- 
las (3.5) and (5.6) in Chapter I). 
They cannot exist in the absence 
of charges and currents. How- 
ever, as wc shall sec just below, 
X Maxwell equations have nonzero 
solutions even in the case of iden- 
tically zero currents and charges 
in the space. Let's study one of 
such solutions. We choose some 
right-oriented rectangular Carte- 
sian system of coordinates and 
take some constant vector k directed along a;-axis (see Fig. 5.1). 
Then we choose another constant vector Aq directed along y-axis 
and consider the following two functions: 




Fig. 5.1 



(5.1) 



A = Aq sin(A; x — cot), 



ip = 0. 



Here k = |k|. Suppose p = and j = 0. Then, substituting (5.1) 
into (4.4) and into Maxwell equations (4.5), we get 



(5.2) 



UJ 

c ' 



It is not difficult to satisfy this condition (5.2). If it is fulfilled, 
then corresponding potentials (5.1) describe plane electromag- 
netic wave, Lu being its frequency and k being its wave-vector, 
which determines the direction of propagation of that plane wave. 
Rewriting (5.1) in a little bit different form 



(5.3) 



A = Aq sin(fc(a; — ct)), 
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we see that the velocity of propagating of plane electromagnetic 
wave is equal to constant c (see (1.5) in Chapter I). 

Now let's substitute (5.1) into (3.4) and calculate electric and 
magnetic fields in electromagnetic wave: 

E = Eq cos{k X - CO t) , Eo = |k|Ao, 

(5.4) H = HoCos(A:a;-a;t), Ho = [k,Ao], 

|Eo| = |Ho| = |k||Ao|. 

Vectors k, Eg, and Hq arc perpendicular to each other, they form 
right triple. Wave (5.4) with such vectors is called plane linear 
polarized electromagnetic wave. Vector Eq is taken for polarization 
vector of this wave. Wave 

E = Eq cos(fca; — ojt) + Hq sin(A;a; — cot), 
H = Ho cos(fc X — oot) — Eq sin(fc x — cot) 

is called circular polarized wave. It is superposition of two linear 
polarized waves. Natural light is also clcctromagnc'tic wave. It 
has no fixed polarization, however it is not circular polarized as 
well. Natural light is a superposition of numerous plane linear 
polarized waves with chaotically distributed polarization vectors. 

§ 6. Emission of electromagnetic waves. 

Plane wave (5.4) is an endless wave filling the whole space. It 

is certainly kind of idealization. Real electromagnetic waves fill 
only some restricted part of the space. Moreover, they are not 
eternal in time: there are sources (radiators) and absorbers of 
electromagnetic fields. Formula (5.4) is an approximate descrip- 
tion of real electromagnetic field in that part of space which is far 
apart from radiators and absorbers. 
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In this section we consider process of generation and radiation 
of electromagnetic waves. Usually radiator is a system of charges 
and currents, which is not static. We describe it by means of 
functions p{r,t) and j(r,t). Let's consider Maxwell equations 
transformed to the form (4.5). These are non-homogeneous differ- 
ential equations. Their solutions are not unique: to each solution 
already found one can add arbitrary solution of corresponding 
homogeneous equations. However, if wo assume p{r,t) and j(r, t) 
to be fast decreasing as r — *• oo and apply similar condition to 
^p{r.t) and A(r, f), then wc restrict substantially the freedom in 
choosing solutions of the equations (4.5). In order to find one 
of such solutions we need fundamental solution of d'Alambert 
operator. This is distribution of the form: 

(6.1) u{r,t) = ^eit)SicH'-\r\'), 

where 9 and S are Heaviside theta-function and Dirac delta- 
function respectively. Function (6.1) satisfies dAlambert equation 
with distribution in right hand side: 

au = 6{t)6{r). 

In physics such objects are called Green functions. Knowing 
fundamental solution (6.1) of d'Alambert operator, now we can 
write solution for the equations (4.5) in form of contractions: 



(6.2) ip = 4:Tru*p, A = — u*}. 

Here * denotes contraction of two distributions, see [1]. Due to 
the properties of this operation from charge conservation law (sec 
formula (5.4) in Chapter I) we derive Lorentzian gauge condition 
(4.4) for scalar and vectorial potentials (6.2). For smooth and 
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sufficiently fast decreasing functions p{r,t) and j(r, potentials 
(6.2) are reduced to the following two integrals: 



Here the quantity r = T(r, f ) is called time delay. It is determined 
by the ratio r = |r — f|/c. Potentials (6.3) are called retarded 

potentials. 

Retarded potentials have transparent physical interpretation. 
Scalar potential (p at the point r at time instant t is a superposi- 
tion of contributions from charges at various points of the space, 
the contribution from the point f being determined not by charge 
density at present time instant t, but at previous time instant 
t — T. Time delay r is exactly equal to the time required for the 
signal spreading with light velocity c from the source point f to 
get to the observation point r. Similar time delay is present in 
formula for vector potential A. 

Note that fundamental solution of d'Alambert equation is not 
unique. For example there is a solution obtained from (6.1) 
by changing t for — r. Such solution corresponds to advanced 
potentials. However, in physics advanced potentials have no 
meaning, since they would break causality principle. 

Let's consider system of charges located in some small domain 
n surrounding the origin. Let R be maximal linear size of 
this domain fl. Using formulas (6.3), we calculate we calculate 
electromagnetic field of the system of charges at the point r which 
is far distant from the domain fi, i.e. |f | < -R <C |r|. Due to these 
inequalities the ratio f/|r| is small vectorial quantity. Therefore 




(6.3) 
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we have the following asymptotic expansions for |r — f | and t — t: 

(r r) 

|r-f| = |r|-^ + ..., 

(6.4) 

f r f + 

t — T = t 1 7—. h . ■ . . 

c |r| c 

The ratio |r|/c in (6.4) determines the time required for electro- 
magnetic signal to get from the domain Q. to the observation point 
|r|. Posterior terms in the series for f — r are estimated by small 
quantity R/c. This is the time of propagation of electromagnetic 
signal within domain f2. 

Denote t' = t — \t\/c and let t — t = t' + 6. For the quantity 9 
we have the estimate \9\ < R/c. Then let's consider the following 
Taylor expansions for p and j : 

p(r,t-r)=p{r,t') + ^^^9 + ..., 

(6-5) 

j(f,f-r)=j(f,0 + ^%^^ + .... 

The condition R <C |r| is not sufficient for the expansions (6.5) to 
be consistent. Use of expansions (6.5) for approximating p{r,t — T) 
and j(f , t — t) is possible only under some additional assumptions 
concerning these functions. Denote by T some specific time for 
which functions p and j within domain f2 change substantially. In 
case when one can specify such time T, the following quantities 
are of the same order, i. e. equally large or equally small: 

dt ' 

(6.6) 
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Now (6.5) can be rewritten as follows: 

p{r,t - t) = p{r,t') +T - + . . . , 

(6.7) 

Correctness of use of expansions (6.7) and (6.5) is provided by 
additional condition R/c<^T. This yields G/T <^1. 

The condition R/c <^ T has simple meaning: the quantity 
u! = 2w/T is a. frequency of radiated electromagnetic waves, while 
A = 2ttc/uj = cT is a wavelength. Hence condition R/c <^ T 
means that wavelength is mach greater than the size of radiator. 

Suppose that both conditions R^cT and R <C |r| are fulfilled. 
Let's calculate retarded vector potential A in (6.3) keeping only 
first term in the expansion (6.5): 

(6.8) A= /jM)d3j; + ___ 

J |r|c 
n 

In order to transform integral in (6.8) let's choose some arbitrary 
constant vector e and consider scalar product (A, e). Having 
defined vector a and function /(f) by the relationships 

a = ^ = grad/, /(f) = (a, f), 

we make calculations analogous to that of (7.5) in Chapter I: 

y"(j,grad/) £r = J diy{fi)£r- j f divider = 

dn on 
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Since e is arbitrary vector, for vector potential A from (6.9) we 
derive the following formula: 

./ ot rc rc 



Here D = D(t') is time derivative of dipole moment D of the 

system of charges at time instant t'. 

In a similar way, keeping only initial terms in the expansions 
(6.4) and (6.5), for scalar potential ip in (6.3) we find 



where Q is total charge enclosed within domain fi. This charge 
does not depend on time since domain O is isolated and we have 
no electric current in outer space. 

Let's compare the expressions under integration in (6.10) and 
(6.11) taking into account (6.6). This comparison yields 

lAI « — ^. 

The estimate R/{cT) <C 1 following from R <^ cT means that 
vectorial potential is calculated with higher accuracy than scalar 
potential. Hence in calculating one should take into account 
higher order terms in expansions (6.4) and (6.5). Then 

Q , fdp{r,t') <r, f) 



|r| J at \y\ c 
(6.12) . ^ 

n 
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Calculating integrals in formula (6.12), we transform it to 
Q (D> r) (D, r) 

Potentials (6.10) and (6.13) are retarded potentials of the 
system of charges in dipole approximation. Dependence of p and j 
on time variable t lead to the dependence of D on t' in them. Let's 
consider asymptotics of of these potentials as r — > oo. Thereby 
we can omit last term in (6.13). Then 

(6.14) ^=^ + ^ + -, A=^ + .... 

|r| Irl^'c |r|c 

Now on the base of formulas (3.4) and (6.14) wc find 

asymptotics of electric and magnetic fields at far distance from 
the system of charges. In calculating rot A and gradt/? we take 
into account that t' = t — \r\/c in argument of D(t') is a quantity 
depending on r. This dependence determines leading terms in 
asymptotics of E and H: 

Vectors E and H (more precisely, leading terms in their asymp- 
totics) are perpendicular to each other and both are perpendicular 
to vector r. This situation is similar to that of plane wave. How- 
ever, in present case we deal with spherical wave being radiated 
from the origin. The magnitude of fields |E| ~ |H| decreases as 
l/|r|, which is slower than for static Coulomb field. Using formula 
(2.5), one can find the density of energy flow for waves (6.15): 

(6.16) S = |^r + .... 
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For modulus of vector S we have |S| ~ l/|rp- This means that 
total flow of energy through the sphere of arbitrarily large sphere 
is nonzero. So we have real radiation of electromagnetic energy. 
The amount of radiated energy is determined by second time 
derivative of dipole moment. Therefore this case is called dipole 
approximation in theory of radiation. 

Exercise 6.1. Applying formula (6.16), find angular distribu- 
tion of the intensity for dipole radiation. Also find total intensity 

of dipole radiation. 

Exercise 6.2. Particle with charge q is moving along circular 
path of radius R with constant velocity v = uR for infinitely long 
time (lo is angular velocity). Calculate retarding potentials and 
find angular distribution for intensity of electromagnetic radiation 
of this particle. Also find total intensity of such cyclotronic radia- 
tion. 

Exercise 6.3. Assume that charge density p is zero, while cur- 
rent density j is given by the following distribution: 

(6.17) j(r,t) = -c[M(i), grad5(r)] 

(compare with (7.16) in Chapter I). Find retarding potentials (6.2) 
for (6.17). Also find angular distribution and total intensity for 
magnetic- dipole radiation induced by current (6.17). 



CHAPTER III 



SPECIAL THEORY OF RELATIVITY 



§ 1. Galileo transformations. 

Classical electrodynamics based on Maxwell equations histor- 
ically was first field theory. It explained all electromagnetic 
phenomena and predicted the existence of electromagnetic waves. 
Later on electromagnetic waves were detected experimentally and 
nowadays they have broad scope of applications in our everyday 
life. However, along with successful development of this theory, 
some difficulties there appeared. It was found that classical elec- 
trodynamics contradicts to relativity principle. This principle in 
its classical form suggested by Galileo and Newton states that 
two Cartesian inertial coordinate systems moving with constant 
velocity with respect to each other are equivalent. All physi- 
cal phenomena in these two systems happen identically and are 
described by the same laws. 

Let's consider two such Cartesian inertial coordinate systems 
{r,t) and {r,i). Suppose that second system moves with velocity 
u relative to first one so that coordinate axes in motion remain 
parallel to their initial positions. The relation of radius-vectors of 
points then can be written in form of the following transforma- 
tions known as Galileo transformations: 

(1.1) t = i, r = f-hui. 
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First relationship (1.1) means that watches in two systems are 
synchronized and tick synchronously. Let r{i) be trajectory of 
some material point in coordinate system (f,f). In first coordi- 
nate system this trajectory is given by vector r(t) = r{i) + ui. 
Differentiating this relationship, due to i = t in (1.1) we get 

„x dr dr 

(1.2) at = at+"' " = " + 

Last relationship in (1.2) is known as classical law of velocity 

addition. Differentiating (1.2) once more, wc find the relation for 
accelerations of material point in these two coordinate systems: 

According to Newton's second law, acceleration of material point 
is determined by force F acting on it and by its mass ma = F. 
Prom (1.3) due to relativity principle we conclude that force F 
is invariant quantity. It doesn't depend on the choice of inertial 
coordinate system. This fact is represented by the relationship 

(1.4) F(f + ut,v + u) =F(f,v). 



Now let's consider charged particle with charge q being at 
rest in coordinate system {r,i). In this coordinate system it 
produces Coulomb electrostatic field. In coordinate system (r, t) 
this particle is moving. Hence it should produce electric field and 
magnetic field as well. This indicate that vectors E and H are 
not invariant under Galileo transformations (1.1). Even if in one 
coordinate system we have pure electric field, in second system we 
should expect the presence of both electric and magnetic fields. 
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Therefore transformation rules for E and H analogous to (1.4) for 
F should be written in the following form: 

E(f + ui, i) = a(E(f , i), H(f , i), u), 
H(f + ui, i) = /3(E(f , i) , H(f , i) , u) . 

Duo to superposition principle, which is fulfilled in both coordi- 
nate systems, functions a and /3 are linear and homogeneous with 
respect to E and H. Therefore (1.5) is rewritten as 



(1.6) 



E(r, t) = ai E(f , i) + a2 H(f , i), 
H(r,<) =/3iE(f,t)+/32H(f,f), 



where ai, a2, /3i, P2 are some linear operators which depend on 
u only. Vectors E and H determine the action of electromagnetic 
field upon charges in form of Lorontz force (see formula (4.4) in 
Chapter I). Substituting (1.6) into that formula and taking into 
account (1.2) and (1.4), we get 

qai E + ga2 H + - [v + u, /3i E]-|- 

(1.7) ^ 

+ ^[v + u, /?2H] = gE+^[v, H]. 
c c 

The relationship (1.7) is an identity with three arbitrary parame- 
ters: V, E, H. Therefore we can equate separately terms bilinear 
with respect to v and E. This yields [v, /3i E] = 0, hence (3i = 0. 
Now let's equate terms bilinear with respect to v and H. This 
yields [v, /?2H] = [v, H]. Hence /32 = 1. And finally we should 
equate terms linear with respect to H and E. This yields the 
following formulas for operators ai and a2- 

a2H=--[u, Hi, ai = 1. 

c 
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Now, if we substitute the above expressions for operators ai, a2, 
01, 02 into formula (1.6), we get the relationships 

(1.8) E = E-i[u, H], H = H. 

The relationships (1.8) should complete Galileo transforma- 
tions (1.1) in electrodynamics. However, as we shall see just 
below, they cannot do this mission in non-contradictory form. 
For this purpose, let's transform Maxwell equations written as 
(1.1) and (1.2) in Chapter II to coordinate system {r,i). For 
partial derivatives due to transformations (1.1) we have 

fe=i 

Now, combining (1.8) and (1.9), we derive 

divH = divH, 

divE = divE+ i(u, rotH), 
c ^ ' 

rotH = rotH 

rotE = rotE-l- -{u, H} - - u divH, 

c c 

an an ~ 

^ = ^ - {u, E} + - [u, {u, H}] - - [u, dil/di]. 



Here by curly brackets we denote commutator of vector fields (see 
[2]). Thereby vector u is treated as constant vector field. 
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When substituting the above expressions into Maxweh equa- 
tions we consider the case of zero charges and currents: p = 0, 
j = 0. This yields the following equations: 

divH = 0, 

divE = --(u, rotH), 
c ^ ' 

- 1 9E 1 , fi, 
rotH= --^--{u, E}+ 

c dt c 

+ l[u, {u, H}]-^[u, dU/di], 

~ 1 ail 

rotE = zr. 

c dt 

Only two of the above four equations coincide with original 
Maxwell equations. Other two equations contain the entries 
of vector u that cannot be eliminated. 

This circumstance that we have foimd is very important. In the 
end of XlX-th century it made a dilemma for physicists. The way 
how this dilemma was resolved had determined in most further 
development of physics in XX-th century. Indeed, one had to 
make the following crucial choice: 

(1) to admit that Maxwell equations are not invariant with 
respect to Galileo transformations, hence they require the 
existence of some marked inertial coordinate system where 
they have standard form given in the very beginning of 
Chapter II; 

(2) or to assume that formulas (1.1) are not correct, hence 
relativity principle claiming equivalence of all inertial co- 
ordinate systems is realized in some different way. 
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First Choice had lead to ether theory. According to this theory, 
marked inertial coordinate system is bound to some hypothetical 
matter, which was called ether. This matter has no mass, no 
color, and no smell. It fills the whole space and does not reveal 
itself otherwise, but as a carrier of electromagnetic interaction. 
Specified properties of ether look quite unusual, this makes ether 
theory too artificial (not natural). As a compromise this theory 
was admitted for a while, but later was refuted by experiments of 
Michaelson and Morley, who tried to measure the Earth velocity 
relative to ether (ether wind). 

Second choice is more crucial. Indeed, refusing formulas (1.1), 
we refuse classical mechanics of Newton in whole. Nevertheless 
the development of science went through this second choice. 



Having refused formulas (1.1), one should replace them by 
something else. This was done by Lorentz. Following Lorentz, 
now we replace Galileo transformations (1.1) by general linear 
transformations relating {r,t) and {r,i): 



In (2.1) we introduced c as a factor for time variables t and t in 
order to equalize measure units. Upon introducing this factor all 
components of matrix 5* appear to be purely numeric quantities 
that do not require measure units at all. It is convenient to denote 
ct by r^ and treat this quantity as additional (fourth) component 
of radius- vector: 
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3 



3 



(2.1) 





(2.2) 



r° = ct. 
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Then two relationships (2.1) can be united into one relationship: 

(2.3) r^^^^^r^ 

fe=0 

In order to have invertible transformation (2.3) one should assume 
that detS ^ 0. Let T = S^^. Then inverse transformation for 

(2.3) is written as follows: 

3 

(2.4) r-=^T^r^ 

fe=0 

By their structure transformation (2.3) and (2.4) coincide with 
transformations of coordinates of four-dimensional vector under 
the change of base. Soon we shall see that such interpretation 
appears to be very fruitful. 

Now the problem of deriving Lorentz transformations can be 
formulated as problem of finding components of matrix S in (2.3). 
The only condition wc should satisfy thereby is the invariancc 
of Maxwell equations with respect to transformations (2.3) upon 
completing them with transformations for /9, j, E and H. 

For the beginning let's consider the case with no currents and 
charges, i. e. the case p = 0, j = 0. Instead of Maxwell equations 
let's study their differential consequences written in form of the 
equations (4.6) in Chapter II: 

(2.5) DE = 0, DH = 0. 

Invariancc of (2.5) under the transformations (2.3) and (2.4) is 
necessary (but possibly not sufficient) condition for invariancc of 
Maxwell equations from which the equations (2.5) were derived. 
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Further we need the following formula for d'Alambert operator 
used in the above equations (2.5): 

i=0 j=0 

Here by we denote components of matrix 
(2.7) g*^ = Qij 



/I 
0-10 

0-10 

Vo -1, 



It is easy to see that inverse matrix gij for (2.7) has the same 
components, i.e. gij = g*-'. 

From (2.3) and (2.4) we derive the following transformation 
rules for first order differential operators: 

k=0 fe=0 

Substituting (2.8) into formula (2.6), we get 

□ = V Vfl^'«— — 

p=0 q=0 

where matrices g'-' and g^'' are related by formula 
(2.9) r = i212Tf^j9''- 

j=0 j=0 
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In terms of inverse matrices Qpq and Qpq this relationship (2.9) can 
be rewritten as follows: 

3 3 

(2.10) 5,,=^^TfT/5p,. 

p=0 q=0 

Theorem 2.1. For any choice of operator coefEcients ol\, a-i, 
(3i, and /32 in formulas (1.6) the invariancc of the form of equations 
(2.5) under the transformations (2.3) and (2.4) is equivalent to 
proportionality of matrices g and g, i. e. 

(2.11) = Xg''- 

Numeric factor A in formula (2.11) is usually chosen to be equal 
to unity: A = 1. In this case from (2.10) and (2.11) we derive 

(2.12) g,,^Y.Y^TfT]gpq. 

p=0 q-O 

In matrix form this relationship (2.12) looks like 

(2.13) T*gT = g. 

Here 5 is a matrix of the form (2.7), while by T* in (2.13) we 
denote transposed matrix T. 

Definition 2.1. Matrix T satisfying the relationship (2.13) is 
called Lorentzian matrix. 

It is easy to check up that the set of Lorentzian matrices form 
a group. This group is usually denoted by 0(1,3). It is called 
matrix Lorentz group. 
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From the relationship (2.13) for Lorcntzian matrix we derive 
the equahty (detT)^ = 1. Hence detT = ±1. Lorentzian matrices 
with unit determinant form the group SO (1,3), it is called special 
matrix Lorentz group. 

If i = j = 0, from (2.12) we obtain the following formula 
relating components of Lorentzian matrix T: 

(2.14) {T°r - {T^f - {Tif - {Tif = 1. 

Inequality |To | ^ 1 is immediate consequence of the relationship 
(2.14). Hence ^ 1 or < -1. Lorcntzian matrix with ^ 1 
is called orthochronous. The set of orthochronous Lorentzian ma- 
trices form orthochronous matrix Lorentz group ©"""(l, 3). Intersec- 
tion S0^(1, 3) = S0(1, 3)nO'''(l, 3) is called special orthochronous 
matrix Lorentz group. 

Exercise 2.1. Prove theorem 2.1 under the assumption that 
transformation (1.6) given by operator coefEcients ol\, a-i, /Ji, and 
/32 is invertible. 

§ 3. Minkowsky space. 

In previous section we have found that each Lorentzian matrix 
from group 0(1,3) determines some transformation (2.1) pre- 
serving the form the equations (2.5). In deriving this fact wc 
introduced notations (2.2) and united space and time into one 
four-dimensional "space-time". Let's denote it by M. Four- 
dimensional space M is basic object in special theory of relativity. 
Its points are called events. The space of events is equipped with 
quadratic form g with signature (1,3). This quadratic form is 
called Minkowsky metric. Thereby incrtial coordinate systems are 
interpreted as Cartesian coordinates for which Minkowsky metric 
has canonical form (2.7). 
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Equivalence principle. All physical laws in any two inertial 

coordinate systems are written in the same form. 

Let's choose some inertial coordinate system. This choice 
determines separation of event space M into geometric space V 
(space of points) and time axis T: 

(3.1) M = T®V. 

Matrix of Minkowsky metric in chosen coordinate system has 
canonic form (2.7). Therefore geometric space V is orthogonal to 
time axis T with respect to Minkowsky metric g. Restriction of 
this metric to V is negative quadratic form. Changing its sign, 
we get positive quadratic form. This is standard Euclidean scalar 
product in V. 

Now let's consider another inertial coordinate system. Like 

(3.1) , it determines second expansion of M into space and time: 

(3.2) M = f®V. 

In general time axes T and T in expansions (3.1) and (3.2) do 

not coincide. Indeed, bases of these two coordinate systems are 
related to each other by formula 

3 

(3.3) ei = Y,Siej, 

3=0 

where S is Lorentzian matrix from (2.3). For base vector eo 
directed along time axis T from (3.3) we derive 

(3.4) eo = ^0° eo + ei + ea + eg . 

In general components Sq, Sq, and Sq in Lorentz matrix S are 
nonzero. Therefore vectors §0 and eo are non-coUinear. Hence 
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Non-coinc;i(lcnce of time axes T for two inertial coordinate 
systems leads to non-coincidence of geometric spaces: V ^ V. 
This fact lead to quite radical conclusion when we interpret 
it physically: observers in two such inertial systems observe two 
different three-dimensional geometric spaces and have two different 
time ticks. However, in our everyday life this difference is very 
small and never reveals. 

Let's calculate how big is the difference in the rate of time ticks 
for two inertial coordinate systems. From (2.4) we get 



Let t — > +00. If Lorentzian matrix T is orthochronous, then 

Tq > and t — > +oo. If matrix T is not orthochronous, then 
t — > +00 we get t ^ —oo. Transformations (2.4) with non- 
orthochronous matrices T invert the direction of time exchanging 
the future and the past. It would be very intriguing to have such 
a feature in theory. However, presently in constructing theory of 
relativity one uses more realistic approach. So we shall assume 
that two physically real inertial coordinate systems can be related 
only by orthochronous Lorentz matrices from 0+(l,3). 

Restriction of the set of admissible Lorentz matrices from 
0(1,3) to 0+(l,3) is due to the presence of additional structure 
in the space of events. It is called polarization. Let's choose some 
physical inertial coordinate system. Minkowsky metric in such 
system is given by matrix of canonical form (2.7). Let's calculate 
scalar square of four-dimensional vector x in Minkowsky metric: 



(3.5) 




fc=i 



(3.6) 



9(x, X) = - {x'f - {x^f - {x^f. 



By value of their scalar square gf(x, x) in Minkowsky metric g 
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vectors of Minkowsky space M are subdivided into three parts: 

(1) tome-like vectors, for which g(x, x) is positive; 

(2) light vectors, for which g{x, x) = 0; 

(3) space-like vectors, for which g{x, x) is negative. 

Coordinates of Hght vectors satisiy the following equation: 

(3.7) {x°f - {x'f - {x^f - {x'f = 0. 

It is easy to see that (3.7) is the equation of cone in four- 
dimensional space (see classification of quadrics in [4]). This cone 
(3.7) is called light cone. 

Time-like vectors fill interior of 
light cone, while space-like vectors 
fill outer space outside this cone. 
Interior of light cone is a union of 
two parts: time-like vectors with 
x^ > are directed to the future, 
others with xq < arc directed to 
the past. Vector directed to the 
future can be continuously trans- 
formed to any other vector directed 
to the future. However, it cannot 
Fig- 3.1 be continuously transformed to a 

vector directed to the past without 
making it space-like vector or zero vector at least once during 
transformation. This means that the set of time-like vectors is 
disjoint union of two connected components. 

Definition 3.1. Geometric structure in Minkowsky space M 

that marks one of two connected components in the set of time- 
like vectors is called polarization. It is used to say that marked 
component points to the future. 
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Let eo, ei, 82, be orthonormal base in Minkowsky metric*. 
In the space M with polarization one can consider only those 
such bases for which unit vector of time axis eo is directed to the 
future. Then transition from one of such bases to another would 
be given by orthochronous matrix from group 0''~(1, 3). 

Definition 3.2. Four-dimensional affine space M equipped 
with metric g of signature (1,3) and equipped with orientation** 
and polarization is called Minkowsky space. 

According to special theory of relativity Minkowsky space, 
which is equipped with orientation and polarization, is proper 
mathematical model for the space of real physical events. Now 
we can give strict mathematical definition of inertial coordinate 
system. 

Definition 3.3. Orthonormal right inertial coordinate system 
is orthonormal right coordinate system in Minkowsky space with 
time base vector directed to the future. 

It is easy to verify that any two inertial coordinate systems as 
defined above are related to each other by Lorentz transforma- 
tion with matrix S from orthochronous Lorentz group SO~^{l, 3). 
Let's choose one of such coordinate systems and consider related 
expansion (3.1). It is clear that eo G T, while linear span of 
spatial vectors ei, e2, defines subspace V. Taking orthonormal 
base ei, 62, ea for the standard of right bases in V, we equip 
this three-dimensional space with orientation. This is concor- 
dance with the fact that geometric space that we observe in our 
everyday life possesses orientation distinguishing left and right. 



* i.e. base for which Minlcowsky metric has the form (2.7). 
** remember that orientation is geometric structure distinguishing left and 
right bases (see [4]). 
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Exercise 3.1. By analogy with definition 3.3 formulate the 
definition of skew-angular inertial coordinate system. 



Galileo transformations are used in mechanics for describing 
physical processes as they are seen by two observers representing 
two inertial coordinate systems. Lorentz transformations, which 
we have derived from the condition of invariance of electrodynam- 
ical equations (2.5), are designed for the same purpose. However, 
this is not immediately clear when looking at formulas (2.3) and 
(2.4). Therefore we shall bring these formulas to the form more 
convenient for studying their physical nature. 

Let's fix two inertial coordinate systems related by Lorentz 
transformation (2.1). First one is related with orthonormal base 
eg, ei, 62, e;j in Minkowsky space and with the expansion (3.1). 
Second is related with the base eo, ei, e2, and with the 
expansion (3.2). If time axes are parallel eo — eo, then Lorentz 
matrix S in (2.3) is reduced to orthogonal matrix O G 50(3) 
relating two right orthonormal bases ei, 62, 63 and §1, §2, §3. It 
has the following blockwise-diagonal shape: 



Thus, in case if T || T two inertial coordinate systems differ only 
in directions of spatial axes. They do not move with respect to 
each other. 

Now let's consider the case T )/( T. Hence eo ^ Bq. Let H be 
linear span of vectors eg and eg. Denote by W intersection of 
subspaces V and V from (3.1) and (3.2): 



§ 4. Kinematics of relative motion. 



(4.1) 



S = 



/I \ 

01 o\ 01 

01 01 ol 

Vo ol ol ol) 



(4.2) 



H = Span(eo,eo) 



w = vr\V. 
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Lemma 4.1. Two-dimensional subspaces H and W in (4.2) are 
perpendicular to each other in Minkowsky metric g. Their inter- 
section is zero: H (IW = {0}, while direct sum of these subspaces 
coincides with the whole Minkowsky space: H ®W = M. 

Proof. Subspacc H is two-dimensional since it is linear span 
of two non-coUinear vectors. Subspaces V and V are three- 
dimensional and V ^ V. Hence their sum V -\-V coincides with 

M , i. e. dim(F + V^) = 4. Applying theorem on the dimension of 
sum and intersection of two subspaces (sec [4]), we get 

dim(T^) = dimy -h dimV - dim(F -hy)=3-h3-4 = 2. 

In order to prove orthogonality of H and W we use orthogonal- 
ity of T and V in the expansion (3.1) and orthogonality of T and 

V in (3.2). Let y be an arbitrary vector in subspace W. Then 
y £V . From V J- T wc get y _L Gq. Analogously from y <E V we 
get y _L eo. Now from orthogonality of y to both vectors eo and 
§0 we derive orthogonality of y to their linear span: y ± H. Since 
y is arbitrary vector in W, we have W -L H. 

Now let's prove that H nW = {0}. Let's consider an arbitrary 
vector X G H OW. From x G H and x G.W due to orthogonality 
of H and W, which is already proved, we get g(x, x) = 0. But 
X £ W C V, while restriction of Minkowsky metric to subspace 

V is negative quadratic form of signature (0,3). Therefore from 
g{x, x) = we derive x = 0. Proposition H DW = {0} is proved. 

From H nW = {0} we conclude that sum of subspaces H and 
W is direct sum and dim{H-\-W) = 2-1-2 = 4. Hence H®W = M. 
Lemma is proved. □ 

Now let's return back to considering pair of inertial coordinate 
systems with bases eo, ei, 62, 63 and Bq, ei, 62, 63. There is the 
expansion (3.4) for vector §0. Let's write it as follows: 



(4.3) 



eo = 5° eo + V. 
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Here v = ei + Sq 62 + Sq 63 e V. Since matrix S is or- 
thochronous and since §0 ^ eo, we have 

(4.4) S° > 1, V ^ 0. 

For any real number a > 1 there exists a number a > such that 
a = cosh(a). Let's apply this observation to Sq in (4.3): 

(4.5) 5g = cosh(a). 

From (4.3), from (4.5), and from orthogonality of vectors eo and 
V in Minkowsky metric we obtain 

1 = g{eo,eo) = {S^ f g{eo,eo) + g{v,v) = cosh^(a) - |vp. 

Using this equality we can find Euclidean length of vector v in 
three-dimensional subspace V: 

(4.6) |v| = sinh(Q;), where a > 0. 

Let's replace vector v by vector of unit length hi = v/|v| and 
rewrite the relationship (4.3) as follows: 

(4.7) eo = cosh(a) eo + sinh(a) hi. 

Due to (4.7) vector hi is linear combination of vectors Gq and eg, 
hence hi G H. But hi gV as well. Therefore hi € VdH. Vectors 
eo and hi are perpendicular to each other, they form orthonormal 
base in two-dimensional subspace H: 

(4.8) (?(eo,eo) = l, 5(hi,hi) = -l. 

From (4.8) we conclude that restriction of Minkowsky metric to 
subspace H is metric with signature (1, 1). 
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Now we need another vector from subspace H. Let's determine 
it by the foUowing relationship: 

(4.9) hi = sinh(a) eo + cosh(a) hi . 

It is easy to check that vectors eo and hi form another orthonor- 
mal base in subspace H. Transition matrix relating these two 
bases has the following form: 

Ccosh(a) sinh(a) \ 
sinh(a) cosh(a) / 

Matrix (4.10) is called the matrix of Lorentzian rotation. 

There is four-dimensional version of matrix (4.10). Indeed, 

vector hi G y is perpendicular to subspace W CV. Therefore we 
have the expansion of subspace F as a direct sum: 



V = Span(hi) e W. 



Let's choose two vectors h2 and hs forming orthonormal base 
in subspace W and complementing hi up to an orthonormal 
right base in V. Then four vectors eo, hi, h2 constitute 
orthonormal right base in M with time vector eo directed to 
the future. Transition matrix relating this base with the base 
eo, hi, h2 ha has the following bio ckwise- diagonal form: 



(4.11) 



( cosh(Q!) sinh(Q!) \ 

sinh(Q:) cosh(a) 

10 

V 1 / 
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Transition from base eo, ei, 62 63 to base eg, hi, h2 ha is given 
by a matrix of the form (4.1). This is because their time vectors 
do coincide. In a similar way transition from base eo, hi, h2 
to base eg, ei, §2 63 is given by a matrix of the same form 
(4.1). Ultimate change of base eo, ei, e2 63 for another base 
§0, ei, §2 63 then can be done in three steps. 

Theorem 4.1. Each Lorentzian matrix S G 50+ (1,3) is a 
product of three matrices S = Si 81,82, one of which Sl is a 
matrix of Lorentzian rotation (4.11), while two others 81 and 82 
are matrices of the form (4.1). 

In order to clarify physical meaning of Lorentz transformations 
let's first consider transformations with matrix 8 of the form 
(4.11). Let ct = r°, , , r"' be coordinates of some vector 
r G M in the base eo, hi, h2, h^. By ci = f°, f ^, f^, we denote 
coordinates of the same vector in the base §0, hi, h2 h^. For 
matrix 8 of the form (4.11) formula (2.3) leads to relationships 

, , , ~ sinh(a) 

t = cosh(a) t -\ — r\ 

c 

= sinh(a) cf + cosh(a) f^, 

r2 = 

Let r^, f^, be coordinates of radius- vector of some point A 
which is at rest in incrtial coordinate system with base ei, 62, §3. 
Then f^, f^, are constants, they do not depend on time t in 
this coordinate system. Upon calculating coordinates of this point 
A in other incrtial coordinate system by means of formulas (4.12) 
its first coordinate appears to be a function of parameter t. 
We use first relationship (4.12) in order to express parameter i 



(4.12) 
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through time variable t in second coordinate system: 

(4.13) ^-^ t^_tanh(a).,_ 

cosh(a) c 

Substituting (4.13) into other three formulas (4.12), we get 

= (t) = c tanh(a) t + const, 

(4.14) = r2(t) = const, 
= r^{t) = const. 

From (4.14) wc sec that in second coordinate system our point A 
is moving with constant velocity u = c tanh(Q;) in the direction of 
first coordinate axis. 

In contrast to parameter a in matrix (4.11), parameter u 
has transparent physical interpretation as magnitude of relative 
velocity of one coordinate system with respect to another. Let's 
express components of matrix (4.11) through u: 

cosh(a) = — , sinh(a) ^ 



Let's substitute these formulas into (4.12). As a result we get 



& 1 
r 



V? 



(4.15) 
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Denote for a while by r and f the following two three- 
dimensional vectors in subspaces V and V: 

r = hi + h2 + hg, 

(4-16) 1-0 , 

Then we define linear map 6 : V ^ V determined by its action 
upon base vectors hi, h2, and h^: 

0(hi)=hi, 0(h2) = h2, 0(h3)=h3. 

This map 6 is orientation preserving isometry, since it maps right 
orthonormal base of subspace V to right orthonormal base in 
subspacc V. Using above notations (4.16) and the map 0, we can 
write formulas (4.15) in vectorial form; 



C2 



(4.17) 

. (eu, r) 
6*11 f + V / 6>u ,^ 

u - (Ou, r) 

0r = , ' + f - \ / eu. 

Here u = uhi is vector of relative velocity of second coordinate 
system with respect to first one. Formulas (4.17) are irrespective 
to the choice of bases in subspaces V and V. Therefore they are 
applicable to Lorcntz transformations with special matrix of the 
form (4.11) and to arbitrary Lorentz transformations with matrix 
5 = 5i S'l ^2 (see theorem 4.1). 
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Very often the sign of map 9 realizing isomorphism, of subspaces 
V and V in formulas (4.17) is omitted: 




(4.18) 




Formulas (4.18) represent "conditionally three-dimensional" un- 
derstanding of Lorentz transformations when vectors r and f 
treated as vectors of the same three-dimensional Euclidean space, 
while t and t are treated as scalar parameters. However, accord- 
ing to modern paradigm four-dimensional Minkowsky space is real 
physical entity, not purely mathematical abstraction convenient 
for shortening formulas (compare (2.3) and (4.17)). 

When writing formulas (4.17) and (4.18) in components we 
should expand vectors r and u in the base of one coordinate sys- 
tem, while vector f is expanded in the base of another coordinate 
system. Thereby the difference in the shape of these two formulas 
completely disappears. 

Exercise 4.1. Using expansions (4.16) for vectors r and f , de- 
rive the following formulas: 



Combining these formulas with (4.15), derive formulas (4.17). 



-1 ^ (^u, f) ^ 
|u| ' 



ha f ^ ha = f - 




eu. 
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§ 5. Relativistic law of velocity addition. 



Classical law of velocity addition was first consequence that we 
obtained from Galileo transformations: 



see formulas (1.2). Replacing Galileo transformations by Lorentz 
transformations, now we should derive new relativistic law of 
velocity addition. 

Suppose that vector-function r{t) describes the motion of a 
point A in inertial coordinate system (f , i) and suppose that this 
coordinate system moves with velocity u with respect to other 
inertial coordinate system {r,t). For passing to coordinate system 
(r,t) we use Lorentz transformation given by formulas (4.18). As 
a result we get two functions 



(5.1) 



V = V + u 




(5.2) 



ut + 



u 



r(i) 




(u, Hi)) 

|U|2 



u. 



Let's calculate first derivatives of functions 



(5.2): 



(5.3a) 
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,,,,, dr "+1^^^" . (u,v) 
(5.3b) = , ' ' + V - V u. 

By V we denote the velocity of the point A in coordinates (f , t): 

i,-, dr 

In a similar way by v we denote the velocity of this point in other 
coordinates {r,t). To calculate v we divide derivatives: 

dr ... f dr\ / f dt 
Substituting (5.3a) and (5.3b) into (5.4), we get formula 

(5.5) 



, (u, v) - (u. v) 

U + ^7 TTT-Vl V ^-r TTT- U 



Formula (5.5) is relativistic law of velocity addition. It is mach 
more complicated than classical law given by formula (5.1). How- 
ever, in the limit of small velocities |u| <C c formula (5.5) reduces 
to formula (5.1). 

Exercise 5.1. Derive relativistic law of velocity addition from 
formula (4.17). Explain why resulting formula differs from (5.5). 

§ 6. World lines and private time. 

Motion of point-size material object in arbitrary inertial coor- 
dinate system {r,t) is described by vector-function r{t), where t 



CopyRight © Sharipov R.A., 1997, 2003. 



92 



CHAPTER III. SPECIAL RELATIVITY. 



is time variable and r is radius- vector of point. Four-dimensional 
radius- vector of this material point has the following components: 

(6.1) rO(t) = ci, r\t), r^{t), r^{t). 

Vector-function with components (6.1) determines parametric line 
in Minkowsky space M, this line is called world line of material 
point. Once world line is given, motion of material point is 
described completely. Let's differentiate four-dimensional radius- 
vector (6.1) with respect to parameter t. As a result we get 
four-dimensional vector tangent to world line: 

(6.2) K = (c,r\f^r3) = {c,v^,v'^,v^). 

Last three components of this vector form velocity vector of 
material point. Velocity of most material objects is not greater 
than light velocity: |v| < c. When applied to vector K in (6.2) 
this means that tangent- vector of world line is time-like vector: 

(6.3) 5(K,K) = c2-|v|2>0. 

Definition 6.1. Smooth curve in Minkowsky space is called 
time-like curve if tangent-vector of this curve is time-like vector 
at each its point. 

World lines for most material objects are time-like curves. 

Exception are world lines of photons (light particles) and world 
lines of other elementary particles with zero mass. For them 
|v| = c, hence we get ^(K, K) = 0. 

World line have no singular points. Indeed, even if (/(K,K) = 
0, tangent vector K in (6.2) is nonzero since K*^ = cj^ 0. 

Let's consider world line of material point of nonzero mass. 
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For this line we have the condition (6.3) fulfilled, hence we can 
introduce natural parameter on this line: 



Integral (6.4) yields invariant parameter for world lines. For any 
two points A and B on a given world line the quantity s{B) — s{A) 
docs not depend on incrtial coordinate system used for calculating 
integral (6.4). This quantity is called interval length of the arc 
AB on world line. 

Theorem 6.1. Straight line segment connecting end points of 

an arc on smooth time-Uke curve is a segment of time-Uke straight 
line. Its interval length is greater than interval length of corre- 
sponding arc. 

Let A and B be two successive events in the "life" of material 
point of nonzero mass. The answer to the question what time in- 
terval separates these two events depend on the choice of inertial 
coordinate system from which we obscirvc; the "life" of this ma- 
terial point. So this answer is relative (not invariant). However, 
there is invariant quantity characterizing time distance between 
two events on world line: 



This quantity r in formula (6.5) is called interval of private time 
on world line. 

Concept of private time determine microlocal concept of time 
in theory of relativity. According to this concept each material 
point lives according to its own watch, and watches of different 
material points are synchronized only in very rough way: they 



(6.4) 




to 
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count time from the past to the future. This rough synchroniza- 
tion is determined by polarization in Minkowsky space. Exact 
synchronization of watches is possible only when material points 
come to immediate touch with each other, i. e. when their world 
lines intersect. However, even after such exact synchronization 
in the point of next meeting watches of tho material points will 
show different times. This difference is due to different "life 
paths" between two meetings. 

Concept of private time is illustrated by so-called twins prob- 
lem, well-known from science fiction. Suppose that one of twins 
goes to far-away travel in interstellar spacecraft, while his brother 
stays on the Earth. Which of them will be older when they meet 
each other on the Earth in the end of space voyage. 

The answer is: that one who stayed on the Earth will be older. 
World lines of twins intersect twice. Both intersections occur on 
the Earth, one before travel and other after travel. Its known that 
Coordinate system associated with the Earth can be taken for 
inertial coordinate system with high degree of accuracy (indeed, 
acceleration due to rotation of the Earth around its axis and due 
to orbital rotation around the Sun is not sensible in our everyday 
life). Therefore world line of twin stayed on the Earth is straight 
line. World line of twin in spacecraft is curved. In the beginning 
of travel he accelerates reaching substantial velocity comparable 
with light velocity in the middle of the path. Then he experiences 
backward acceleration in order to brake before reaching target 
point of his travel. Then he accelerates and brakes again in his 
back way to the Earth. According to theorem 6.1 interval length 
of curved world line connecting two events is shorter than interval 
length of straight world line connecting the same two events. 
Hence twin stayed on the Earth will be older. 

Exercise 6.1. Remember proof of the fact that the length of 
curved Une connecting two points A and B in Euclidean space is 
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greater than the length of straight hue segment AB. By analogy 
to this proof find the proof for theorem 6. 1 . 

§ 7. Dynamics of material point. 

Motion of material point in theory of relativity is described by 
its world line in Minkowsky space. Let's choose natural parameter 
on world line and consider four-dimensional tangent vector 

dr(s) 



(7.1) 



u(.) 



ds 



where r(s) is four- dimensional radius vector of events on world 
line. Vector u is (7.1) is called vector of A-velocity. It is time-like 
vector and it is unit vector in Minkowsky metric: ^(u, u) = 1. 
Upon choosing some inertial coordinate system we can write 
components of 4-velocity vector explicitly: 



(7.2) 



Here v 



u 



are components of three-dimensional velocity 



vector V. Note that components u 



of 4-velocity vector 



are absolute numbers (without measure unit). It is easy to see 
from (7.2). Upon multiplying u by scalar mc with the measure 
unit of momentum we get vector of A-momentum 



(7.3) 



1 





c 















for material point with mass m. Vector p plays important role 
in physics since there is fundamental law of nature: the law of 
conservation of A-momentum. 
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Momentum conservation law. Vector of A-momentum of 
material point which do not experience external action remains 
unchanged. 

According to the law just stated, for particle that do not 
experience external action we have p = const. Hence u = const. 
Integrating the equation (7.1), for r(s) we derive 



r(s) = ro + us. 



Conclusion: in the absence of external action material point moves 
uniformly along straight line. 

External actions causing change of 4-momentum of material 
point are subdivided into two categories: 

(1) continuous; 

(2) discrete. 

Continuous actions are applied to material particle by external 
fields. They cause world line to bend making it curved line. In 
this case p 7^ const. Derivative of 4-momentum with respect to 
natural parameter s is called vector of A-force: 

(«) f^FM^ 

Vector of 4-force in (7.4) is quantitative characteristic of the 
action of external fields upon material particle. It is determined 
by parameters of particle itself and by parameters of external 
fields at ciirrent position of particle as well. We know that 
vector of 4- velocity u is unit vector, therefore 5(p,p) = m^c^. 
Differentiating this relationship with respect to s and taking into 
account that components of matrix (2.7) are constant, we find 



(7.5) 



5(u,F) = 0. 
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The relationship (7.5) means that vector of 4-force is perpendicu- 
lar to vector of 4-velocity in Minkowsky metric, i. e. force vector 
is perpendicular to world line of particle. 

Choosing some inertial coordinate system, we can replace nat- 
ural parameter s in (7.5) by time variable t of this coordinate 
system. Then, taking into account (7.3), from (7.4) we derive 



(7.6) = - |v|2 F\ where i = 1,2,3. 

Now, if we denote by f three-dimensional vector with compo- 
nents = \/ — |vp F"^ , then for three-dimensional vector of 
momentum from (7.6) we obtain differential equation 

The equation (7.7) is treated as relativistic analog of Newton's 
second law. Instciad of c;lassical formiila p = rnv relating momen- 
tum and velocity vectors here we have the following relationship: 

(7.8) p = 



In order to write (7.8) in classical form we introduce the quantity 
(7.9) m„ = ™ 



1 - 



Constant m is called mass at rest, while m„ in (7.8) is called 
dynamic mass of moving particle. Now p = v, and Newton's 
second law is written as follows: 



(7.10) 



{my = f . 
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Formulas (7.9) and (7.10) are the very ones which are in mind 
when one says that mass in theory of relativity depends on 
velocity. It seems to me that such terminology is not so good. 
In what follows we shall mostly use four-dimensional invariant 
equation (7.4) and, saying mass, we shall imply mass at rest. 

Discrete external actions appear 




in those situations when 4-momcn- 
tum of material particle changes 
abruptly in jump-like manner. Such 
situation arise in particle collisions, 
particle confluence, and particle de- 
cay. Collision of particles corre- 
spond to that point in Minkowsky 
space where world lines of two or 



several particles come together. Af- 
Fig. 7.1 ter collision particles can simply fly 

out from that point. But if these 
are molecules of ingredients in chemical reaction, then after col- 
lision we would have new molecules of reaction products. In a 
similar way in collisions of atomic nuclei nuclear reactions occur. 

Let's consider simultaneous collision of k particles. Denote by 
Pi, . . . ,pk their 4-momcnta just before the collision. Suppose that 
as a result of collision he have n new particles created from initial 
ones. Denote by pi , . . . , p„ 4-momenta of outgoing particles just 
after the collision. If fc = 1 this is particle decay process, while if 
n = 1 we have particle confluence into one composite particle. 

Momentum conservation law. Total A-momentum of ingo- 
ing particles before collision is equal to total A-momentum of out- 
going particles after collision: 

k n 

(7.11) EP' = Ep- 

i=l i=l 
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As an example we consider process of frontal collision of two 
identical particles of mass m leading to creation of one particle of 
mass M. Suppose that velocities of initial particles are equal by 
magnitude but opposite to each other: 



Pi = 



P2 = 



1 - 



—V" 



For 4-momentum of resulting particle we have 



Pi 



M 



1 



Applying momentum conservation law (7.11) to this situation, we 
get w = and additionally we obtain 



(7.12) 



M 



1 



From (7.12) we see that mass at rest of resulting composite 
particle is greater than sum of rest masses of its components: 
M > m + m. Conclusion: the low of mass conservation is fulfilled 
approximately only in the limit of small velocities |v| <C c. 

Let's multiply zeroth component of 4-momentum of material 
particle by c. Resulting quantity has the measure unit of energy. 
Let's denote this quantity by E: 
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The quantity (7.13) is called kinetic energy of moving particle. 
Writing relationship (7.11) for zeroth components of colliding 
particles, we get energy conservation law: 

k n 

(7.14) j2^^ = J2^i- 

i=l i=l 

Thus, 4-momcntum conservation law for collision includes both 
energy conservation law (7.14) and the law of conservation for 
three-dimensional momentum. 

Note that for zero velocity v = the above quantity (7.13) 
does not vanish, but takes nonzero value 

(7.15) E = mc^. 

This quantity is known as rest energy of material particle. For- 
mula (7.15) is well-known. It reflects very important fact absent 
in classical physics: the ability of energy to mass and mass to 
energy conversion. In practice conversion of energy to mass is 
realized in particle confluence (sec M > m + m in formula (7.12)). 
Converse phenomenon of particle decay yields mass defect (mass 
decrease). Lost mass is realized in additional amount of kinetic 
energy of outgoing particles. Total conversion of mass to en- 
ergy is also possible. This happens in process of annihilation, 
when elementary particle meets corresponding antiparticle. Large 
amount of energy released in annihilation is scattered in form of 
short-wave electromagnetic radiation. 

§ 8. Four-dimensional form of Maxwell equations. 

Starting from electromagnetic equations DE = and DH = 
in previous sections we have constructed and described Lorentz 
transformations preserving form of these equations. We also have 
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given geometric and physical interpretation of Lorentz transfor- 
mations and even have described dynamics of material points on 
the base of new relativistic notion of space and time. Now time 
has come to remember that equations DE = and DH = 
are not primary equations of electrodynamics, they were derived 
from Maxwell equations. To have complete picture we should we 
should write Maxwell equations in four- dimensional form. Let's 
begin with second pair of these equations containing charges and 
currents (see equations (1.2) in Chapter II. Let's modify them: 



1 aE „ An . X. . 

- — rot H = J, - div E = -Awp. 

c at c 



Then rewrite these equations in components using Levi-Civita 
symbol for to express rotor (see [3]): 



(8.1) 



9=1 fe=l 



^ dE'i 

> — — = -47rp. 

9=1 



Here we used notation r° = ct associating time variable with 
zeroth component of radius-vector in Minkowsky space. 

Using Levi-Civita symbol and components of vector H, we can 
construct skew-symmetric 3x3 matrix with elements 



(8.2) 



3 
k=l 



102 



CHAPTER III. SPECIAL RELATIVITY. 



Due to (8.2) we can easily write explicit form of matrix F: 

-H^ 
ifi 



(8.3) 




Let's complement the above matrix (8.3) with one additional line 
and one additional column: 



(8.4) 






-E^ 


-E^ 


-E^ \ 



















\E^ -H^ 



/ 



Additional line and additional column in (8.4) arc indexed by 
zero, i. e. indices p and q run over integer numbers from to 3. 
In addition, we complement three-dimensional vector of current 
density with one more component 



(8.5) 



pc. 



By means of (8.4) and (8.5) we can rewrite Max;well equations 
(8.1) in very concise four-dimensional form: 



i.6) 



q=0 



Now let's consider first pair of Maxwell equations (see equations 
(1.1) in Chapter II). In coordinates they are written as 



(8.7) 



dHP ^ A dE'' 

9=1 fe=l 



dri 



E 

9=1 



dHi 
dri 
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The structure of the equations (8.7) is quite similar to that of 
(8.1). However, their right hand sides are zero and we see shght 
difference in signs. Main difference is that components of vectors 
E and H have exchanged their places. To exchange components 
of vectors E and H in matrix (8.4) we need four- dimensional 
analog of Levi-Civita symbol: 



^pqks 



^pqks 



0, if among p, q, k, s there 
are at least two equal num- 
bers; 

1, ii {pqk s) is even permuta- 
tion of numbers (0123); 

— 1, if (pqk s) is odd permuta- 
tion of numbers (0123). 



Let's define matrix G by the following formula for its components: 



3 3 3 3 
;c=0 s=0 rn=0 n=0 



-^'^ 9km 9sn ^ 



Here g is matrix (2.7) determining Minkowsky metric. Matrix G 
with components (8.8) can be expressed in explicit form: 



(8.9) 



I 

ij2 





-E^ 



E^ 



\ 



E^ -E^ / 



The structure of matrix (8.9) enable us to write remaining two 
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Maxwell equations (8.7) in concise four- dimensional form: 

(-0) t^-. 

Usage of both matrices F and G in theory is assumed to be too 
excessive. For this reason equations (8.10) are written as 

(8-11) EEE-^^"l^-o- 

g=0 fe=0 s=0 

Matrix Fks is obtained from F"*" by means of standard index 
lowering procedure using matrix (2.7): 

3 3 

(8.12) Fks=J2Yl 9km 9sn i^"". 

m=0 n=0 

Four-dimensional form of Maxwell equations (8.6) and (8.11) 
gives a hint for proper geometric interpretation of these equations. 
Matrix (8.4) defines tensor of the type (2, 0) in Minkowsky space. 
This tensor is called tensor of electrowMgnetic field. Tensorial 
interpretation of matrix (8.4) immediately yields transformation 
rules, which were lacking so far: 

3 3 

(8.13) F^"'= ^^^^^S^F™". 

m—O n—0 



These relationships (8.13) determine transformation rules for com- 
ponents of vectors E and H. Before now we express these rules in 
undetermined form by the relationships (1.6). For special Lorentz 
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matrices (4.11) vectors of electric and magnetic fields E and H in 
two inertial coordinate systems are related as follows: 



E'=E\ E^ 



E^ + 


c 


E^ = 






2 ' 










^E' 

c 


H^ = 


H^ + -E^ 




2 ' 




1— ' 





According to theorem 4.1, general Lorcntz matrix is a product 
of special Lorentz matrix of the form (4.11) and two matrices of 
spatial rotation in tree-dimensional space. The latter ones can be 
excluded if one writes Lorentz transformation in "conditionally 
three-dimensional" vectorial form: 

E-il^u-l[u,H] 



(8.14) 



1 - 



(u, H) 



H-^u+i[u, E] 



c 



From (8.13) we derive the following rule for transforming co- 
variant components of the tensor of electromagnetic field: 

3 3 

(8.15) F^,= Y,T.^rT^^rnn. 

m=0 n=0 
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This relationship (8.15) provides invariance of the form of Maxwell 
equations (8.11) under Lorentz transformation (2.3). In order to 
verify this fact it is sufficient to apply relationships (2.8) for 
transforming derivatives and then remember well-known property 
of four- dimensional Levi-Civita symbol ePi^^; 

3 3 3 3 

(8.16) X! ^« Td e'"'""^ = det T s^'^^' . 

a=0 6=0 c=0 d=0 

The condition of invariance of Maxwell equations (8.6) with re- 
spect to Lorentz transformations leads to the following transfor- 
mation rule for components of four-dimensional current density: 

3 

(8.17) f=Y.^^rn~r- 

m— 

In (8.17) it is easy to recognize the transformation rule for com- 
ponents of four-dimensional vector. In case of special Lorentz 
matrix of the form (4.11), taking into account (8.5), one can write 
the above relationship (8.17) as follows: 



(8.18) 




Remember that here u = ctanh(a) is a magnitude of relative 
velocity of one inertial coordinate system with respect to another. 
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In vectorial form relationships (8.18) are written as 



w 

uj) 

|uP " , V <u,j) 

In such form they give transformation rule for charge density p 
and three-dimensional current density j under Lorentz transfor- 
mations with arbitrary Lorentz matrix. 

Exercise 8.1. Prove the relationship (8.16), assuming T to be 

an arbitrary 4x4 matrix. 

Exercise 8.2. Using (2.12), derive the relationship (8.15) from 
(8.12) and (8.13). 

Exercise 8.3. Using (8.15), (8.16) and (2.8), transform Max- 
well equations (8.11) from one inertial coordinate system to an- 
other. Verify that the form of these equations is invariant. 

Exercise 8.4. Using (8.13), (8.17) and (2.8), transform Max- 
well equations (8.6) from one inertial coordinate system to another. 
Verify that the form of these equations is invariant. 

§ 9. Four-dimensional vector-potential. 

Due to special structure of Maxwell equations one can intro- 
duce vector-potential A and scalar potential This was done in 



P + 



(8.19) 
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§ 3 of Chapter II. Here are formulas for components of E and H: 

dffi 1 dAP 



EP = - 



drP c dt ' 



(9.1) 3 3 

jjp _ ^pqk 

g=l fc=l 



^ dri ' 



(see formulas (3.4) in Chapter II). Denote — and consider 
four-dimensional vector A with components A^, A^, A^. This 
is four- dimensional vector-potential of electromagnetic field. By 
lowering index procedure we get covector A: 

3 

(9.2) Ap = Y,9p,A'i. 

q=0 

Taking into account relationships (2.7) for components of matrix 
Qpq, from formula (9.2) we derive 

Aq = Ai = -A\ 

A2 = -A^ As = -A\ 



Moreover, let's write explicitly covariant components for the ten- 
sor of electromagnetic field: 






E^ 


E^ 


E^ \ 


-E^ 







-E^ 









-E^ 






/ 
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Due to (9.3) and (9.4) first relationship (9.1) can be written as 
(9.5) Fo, 



^09 Q^O dj.1 ' 



In order to calculate other components of tensor Fpq let's ap- 
ply (8.2) and second relationship (9.1). Thereby let's take into 
account that Fpg = F^i and Ap = -A^ for p,q= 1, 2, 3: 



(9.6) Fpg = -Y: epg, if'^ = E E E ^-'^ 



dr" 

k=l k=l m=l n=l 



Further transformation of (9.6) is based on one of the well-known 
contraction identities for Levi-Civita symbol: 



(9.7) 5]£p,fce'=""-5;"d--<5;"<5;. 
fe=i 

Applying (9.7) to (9.6), we get 

(9-8) = E E(^p"^," - ^r^.") ^ = ^-'^- 

m—l n—1 

Combining (9.8) and (9.5), we obtain the following formula for all 
covariant components of the tensor of electromagnetic field: 

(99) F 

In essential, formula (9.9) is four-dimensional form of the rela- 
tionships (9.1). It unites these two relationships into one. 

Remember that vectorial and scalar potentials of electromag- 
netic field are not unique. They are determined up to a gauge 
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transformation (see formula (4.1) in Chapter II). This uncertainty 
could be included into transformation rule for component of four- 
dimensional potential A. However, if we assert that A", A^, A^, 
are transformed as components of four-dimensional vector 

3 

(9.10) AP = Y^SPAi, 

g=0 

and Ai, A2, A3 arc obtained from them by index lowering 
procedure (9.2), then we find that quantities Fpq defined by 
formula (9.9) are transformed exactly by formula (8.15), as they 
actually should. 

From (9.10) one can easily derive explicit transformation 
formulas for scalar potential ip and for components of three- 
dimensional vector-potential A. For special Lorentz transforma- 
tions with matrix (4.11) they are written as follows: 



-^1 

A' 



(9.11) 



u 
c 












A\ 





Note that one can rederive transformation rules for components 
of electric and magnetic fields (see §8 above). However, we shall 
not do it now. 

In case of Lorentz transformations with arbitrary Lorentz ma- 
trix the relationships (9.11) should be written in vectorial form: 
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Ill 



(9.12b) 



A = 




u. 



Theorem 9.1. Each skew-symmetric tensor field F of type 
(0, 2) in four-dimensional space satisfying differential equations 
(8.11) is determined by some covector Geld A according to the 
above formula (9.9). 

Proof. Each skew-symmetric tensor field F of type (0, 2) 
in four-dimensional space can be identified with pair of three- 
dimensional vector fields E and H depending on additional pa- 
rameter r° = ct. In order to do this one should use (9.4). Then 
equations (8.11) are written as Maxwell equations for E and H: 



Further construction of covector field A is based on considerations 
from §3 of Chapter II, where three-dimensional vector-potential 
and scalar potential were introduced. Then we denote = 
and thus convert three-dimensional vector-potential into four- 
dimensional vector. And the last step is index lowering procedure 
given by formula (9.2). □ 

Choice of vector field A in formula (9.9), as we noted above, 
has gauge uncertainty. In four-dimensional formalism this fact is 
represented by gauge transformations 



where tp — is some arbitrary scalar field. Formula (9.13) is four- 
dimensional version of gauge transformations (4.1) considered 



divH = 



1 dU 

rotE = — — . 
c at 



(9.13) 
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in Chapter II. It is easy to verify that gauge transformations 
(9.13) do not break transformation rules (9.10) for contravariant 
components of vector A. 

Exercise 9.1. Prove theorem 9.1 immediately in four-dimen- 
sional form without passing back to three-dimensional statements 
and constructions. 

§ 10. The law of charge conservation. 

Earlier wc have noted that charge conservation law can be 
derived from Maxwell equations (see § 1 in Chapter II). To prove 
this fact in four- dimensional formalism is even easier. Let's 
differentiate the relationship (8.6) with respect to and add one 
more summation with respect to index p: 

^ ■ drPdri c ^ dr^ 

p=0 g=0 p=0 

Double differentiation in (10.1) is symmetric operation, while 
tensor of electromagnetic field F^*, to which it is applied, is 
skew-symmetric. Therefore the expression under summation in 
left hand side of (10.1) is skew-symmetric with respect to indices 
p and q. This leads to vanishing of double sum in left hand side 
of formula (10.1). Hence we obtain 

The equality (10.2) is four-dimensional form of charge conserva- 
tion law. If we remember that = cp and = ct, we see that 
this equality coincides with (5.4) in Chapter I. 

Conservation laws for scalar quantities (those like electric char- 
ge) in theory of relativity are expressed by equations analogous 
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to (10.2) in form of vanishing of four-dimensional divergencies for 
corresponding four- dimensional currents. For vectorial quantities 
corresponding current densities are tensors. Thus the law of con- 
servation of 4-momentum for fields is represented by the equation 

p=0 

Tensor T'^p in (10.3) playing the role of current density of 4- mo- 
mentum is called energy-momentum tensor. 

Theorem 10.1. For any vector field j in n-dimensional space 
(n^2) if its divergency is zero 

p=i 

then there is skew-symmetric tensor field ip of type (2, 0) such that 



(10.5) f = Y^ 



dr" ' 

g=l 



Proof. Choosing some Cartesian coordinate system, we shall 
construct matrix ip^'' of the following special form: 











(10.6) 



... ipn-ln 

\— ... -^"-1" / 
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Matrix (10.6) is skow-symmetric, it has (n — 1) independent com- 
ponents. From (10.5) for these components we derive 



— , where fc = 1, . . . , n — 1, 



(10.7) „_ 



/s=l 



Let's define functions V' in (10.7) by the following integrals: 



(10-8) 



+ ^ / r{r\... ,y,... ,r"-\0)dy. 
n-1 J 



It is easy to verify that functions (10.8) satisfy first series of 
differential equations (10.7). Under the condition (10.4) they 
satisfy last equation (10.7) as well. Thus, theorem is proved. □ 

Theorem 10.1 can be easily generalized for arbitrary tensorial 
currents. Its prove thereby remains the same in most. 

Theorem 10.2. For any tensorial field T of type (m, s) in the 
space of dimension 2 if its divergency is zero 
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then there is tensorial field ip of type {m + l,s) skew-symmetric 
in last pair of upper indices and such that 



Exercise 10.1. Verify that the equation (10.4) provides last 
equation (10.7) to be fulfilled for the functions (10.8). 

Exercise 10.2. Clarify the relation of theorem 10.1 and theo- 
rem on vortex field in case of dimension n = 3. 



In previous three sections we have managed to write in four- 
dimensional form all Maxwell equations, charge conservation law, 
and the relation of E, H and their potentials. The relationships 
(8.6), (8.11), (9.9), (9.13), (10.2), which were obtained there, pre- 
serve their shape when we transfer from one rectangular Cartesian 
coordinate system to another. Such transitions arc interpreted as 
Lorentz transformations, they are given by Lorentz matrices. 
However, all these relationships (8.6), (8.11), (9.9), (9.13), (10.2) 
possess transparent tensorial interpretation. Therefore they can 
be transformed to any skew-angular Cartesian coordinate system 
as well. Thereby we would have minor differences: the shape 
of matrix g would be different and instead of gPs*^* in (8.11) we 
would require volume tensor with components 




Pm+l = l 



§11. Note on skew-angular and 
curvilinear coordinates. 



(11.1) 
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Matrix gpg in skew-angular coordinate system is not given by 
formula (2.7), here it is arbitrary symmetric matrix determining 
quadratic form of signature (1,3). Therefore differential equations 
□E = and DH = 0, which we are started from, have not their 
initial form. They are written as □F*'^ = 0, where d'Alambert 
operator is given by formula (2.6) with non-diagonal matrix g^^ . 

In arbitrary skew-angular coordinate system none of axes 
should have time-like direction. Therefore none of them can 
be interpreted as time axis. Three-dimensional form of electrody- 
namics equations, even if we could write them, would not have 
proper physical interpretation in such coordinate system. In par- 
ticular, interpretation of components of tensor F^'^ as components 
of electric and magnetic fields in formula (8.4) would not be 
physically meaningful. 

Tensorial form of four-dimensional electrodynamics equations 
enables us to make one more step toward increasing arbitrariness 
in the choice of coordinate system: we can use not only skew- 
angular, but curvilinear coordinates as well. To make this step we 
need to replace partial derivatives by covariant derivatives: 



(11.2) 




(see [3] for more details). Connection components required for 
passing to covariant derivatives (11.2) are determined by compo- 
nents of metric tensor. The latter ones in curvilinear coordinate 
system do actually depend on r°, r^, r^, r^: 



(11.3) 




dgis 



Qj.s 



) 



No we give list of all basic equations, which we derived above, in 
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covariant form. Maxwell equations are written as follows: 
(11.4) 

3 3 3 
g=0 fe=0 s=0 

Here components of volume tensor ojPi''^ are given by formula 
(11.1). Tensor of electromagnetic field is expressed through four- 
dimensional vector-potential by formula 

(11.5) Fpg = \/pAq - S/qAp, 

while gauge uncertainty in the choice of vector-potential itself is 
described by the relationship 

(11.6) Ak^Ak + Vki^, 

where tp is arbitrary scalar field. Charge conservation law in 
curvilinear coordinates is written as 

3 

(11.7) Y.Vpf = 0. 

p=0 

Instead of formula (2.6) for D'Alambert operator here we have 

3 3 

(11.8) □ = ^^ff*-'V.V,. 

i=0 j=0 

Dynamics of material point of nonzero mass m 7^ is described 
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by ordinary differential equations of Newtonian type: 

F 

(11.9) r = u, VsU= — . 

mc 

Here dot means standard differentiation with respect to natural 
parameter s on world line, while Vg is covariant derivative with 

respect to the same parameter. 

Exercise 11.1. Using symmetry of Christoffel symbols (11-3) 
with respect to lower pair of indices i and j, show that the re- 
lationship (11.5) can he brought to the form (9.9) in curvilinear 
coordinate system as well. 



CHAPTER IV 



LAGRANGIAN FORMALISM 
IN THEORY OF RELATIVITY 



§ 1. Principle of minimal action 
for particles and fields. 

Dynamics of material points in theory of relativity is described 
by their world lines. These are time-like lines in Minkowsky space. 

Let's consider some world line corrc;- 
sponding to real motion of some particle 
under the action of external fields. Let's 
fix two points A and B on this world 
line not too far from each other. Then 
consider small deformation of world line 
in the range between these two points 
A and B. Suppose that we have some 
coordinate system in Minkowsky space 
pj^g J ^ (either Cartesian, or curvilinear, no mat- 

ter). Then our world line is given in 
parametric form by four functions 

(LI) r°(s), r'{s), r^s), r^s), 

where s is natural parameter. Then deformed curve can be given 
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by the following four functions: 

(1.2) f\s) = r\s) + h\e, s), i = 0, . . . , 3. 

Here s is original natural parameter on initial non-deformed world 
line (1.1), while h^{e,s) are smooth functions which are nonzero 
only within the range between points A and B. Note that 
functions /i*(e, s) in (1.2) depend on additional parameter e which 
is assumed to be small. Moreover, we shall assume that 

(1.3) h\e,s) ^0 as e ^ 0. 

Thus, in (1.2) we have whole family of deformed lines. This family 
of lines is called variation of world line (1.1). Due to (1.3) we 
have the following Taylor expansion for h^{e, s): 

(1.4) h'{e,s)=eh\s) + ... . 

Under the change of one curvilinear coordinate system for an- 
other quantities h^{s) are transformed as components of four- 
dimcnsional vector. This vector is called vector of variation of 
world line, while quantities 

(1.5) 5r'{s)=eh\s) 

are called variations of point coordinates. It is clear that they 
also are transformed as components of four-dimensional vector. 
Due to formulas (1.4) and (1.5) parametric equations of deformed 
curves (1.2) are written as follows: 

(1.6) f\s) = r\s) + 5r\s) + . . . . 

By this formula we emphasize that terms other than linear with 
respect to small parameter e are of no importance. 
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By varying functions /i*(e, s) in (1.2) and by varying parameter 
e in them we can surround segment of initial world line by a 
swarm of its variations. Generally speaking, these variations do 
not describe real dynamics of points. However, they are used 
in statement of minimal action principle. Within framework of 
Lagrangian formalism functional of action S is usually introduced, 
this is a map that to each line connecting two points A and B put 
into correspondence some real number S. 

Principle of minimal action for particles. World line con- 
necting two points A and B describes real dynamics of material 
point if and only if action functional S reaches local minimum on 
it among other lines being its small variations. 

Action functional S producing number by each line should de- 
pend only on that line (as geometric set of points in M), but 
it should not depend on coordinate system (r'^ , , r'^ , r^) in M. 
By tradition this condition is called Lorentz invariance, though 
changes of one curvilinear coordinate system by another form 
much broader class of transformations than Lorentz transforma- 
tions relating two rectangular Cartesian coordinate systems in 
Minkowsky space. 

Action functional in most cases is integral. For single point of 
mass m in electromagnetic field with potential A it is written as 



Here q is electric charge of particle, while u = u(s) is vector of 
its 4- velocity (unit tangent vector of world line). First integral 

in (1.7) yields action for free particle (in the absence of exter- 
nal fields), second integral describes interaction of particle with 
electromagnetic field. 



(1.7) 
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If we consider system of A'' particles, then we should write 
integral (1.7) for each of them and we should add all these 
integrals. And finally, in order to get the action functional for 
total system of field and particles we should add integral of action 
for electromagnetic field itself: 

S2(i) S2(«) 

ds I - 



5 = ^ -rui c j ds-^ J g{A, u) 

^2 3 3 



16; 

Vi 



p=0 q=0 




Last integral in (1.8) deserves special consideration. This is 
four-dimensional volume integral over the domain enclosed be- 
future tween two three-dimensional hyper- 

surfaces Vi and V2. Hypersurfaces 
Vi and V2 are space-like, i. e. their 
normal vectors are time-like vectors. 
These hypersurfaces determine the 
fissure between the future and the 
past, and over this fissure we inte- 
grate in (1.8). Thereby change of 
Fig. 1.2 field functions (here these arc com- 

ponents of vector-potential A) when 
passing from V\ to V2 reflects evolution of electromagnetic fleld 
from the past to the future. 

Electromagnetic field is described by field functions. Therefore 
variation of field is defined in other way than that of particles. 
Suppose that f7 is some restricted four-dimensional domain en- 
closed between hypersurfaces V\ and V2. Let's consider four 
smooth functions /i*(e, r) = /i*(£, r°, r^, r^, r^) being identically 
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zero outside the domain CI and vanishing for £ = 0. Let's define 

(1.9) A\r) = A'{r) + h\s,r) 

and consider Taylor expansion of at the point s = 0: 

(1.10) h'{e,r)=eh'{r) + ... . 

The following fnnctions determined by linear terms in the above 
Taylor expansions (1.10) 

(1.11) SA\r)=eh\r) 

are called variations of field functions for electromagnetic field. 
Deformation of vector-potential (1.9) now can be written as 

(1.12) A\r)=A\r) + 6A\r) + .... 

Principle of minimal action for fields. Field functions de- 
termine actual conRguration of physical Gelds if and only if they 
realize local minimum of action functional S in class of all varia- 
tions with restricted support*. 

The condition of minimum of action for actual field configura- 
tion and for actual world lines of particles, as a rule, is not used. 
In order to derive dynamical equations for fields and particles it is 
sufficient to have extremum condition (no matter minimum, max- 
imum, or saddle point). For this reason minimal action principle 
often is stated as principle of extremal action. 

Exercise 1.1. Verify that h^{s) in (1.4) are transformed as 
components of vector under the change of coordinates. 



* Variations with restricted support are those which are identically zero 
outside some restricted domain Q. 
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Exercise 1.2. Prove that under gauge transformations (11.6) 
from Chapter III action functional (1.8) is transformed as follows: 



N 

(1.13) 5 -> 5 - ^ (I Hr{s,m - I V(r(si(z)))) . 



Explain why terms added to action functional in (1.13) are not 
sensitive to variation of word lines (1.2). 

§ 2. Motion of particle in electromagnetic field. 

In order to find world line of relativistic particle in external 
electromagnetic field we shall apply particle version of extremal 
action principle to functional (1.8). Let's choose one of N par- 
ticles in (1.8) and consider deformation (1.6) of its world line. 
When we substitute deformed world line into (1.8) in place of 
initial non-deformed one the value of last integral remains un- 
changed. Thereby in first term containing sum of integrals only 
one summand changes its value, that one which represent the par- 
ticle we have chosen among others. Therefore writing extremity 
condition for (1.8) we can use action functional in form of (1.7). 
The value of (1.7) for deformed world line is calculated as follows: 

(2.1) 5def =-mc J V5(K,K) ds - ^ j g{A,V.)ds. 

Sl Si 

Formula (2.1) visually differs from formula (1.7) because s is 
natural parameter on initial world line, but it is not natural 
parameter on deformed line. Here tangent vector 

(2.2) K(a) = ^=u(s)+s^ + ... 

ds ds 
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is not unit vector. Therefore first integral (1.7) is rewritten as 
length integral (see (6.4) in Chapter III). In second integral (1.7) 
unit tangent vector is replaced by vector K. 

Let's write in coordinate form both expressions which are under 
integration in (2.1) taking into account that we deal with general 
curvilinear coordinate system in Minkowsky space: 



Let's substitute (2.2) into (2.3) and take into account (1.2) and 
the expansion (1.4). As a result for the expressions (2.3) we get 
the following power expansions with respect to small parameter e: 



3 3 



(2.3) 




3 



i=0 



Vai^, u) 



£ 




dh'js) 
ds 



+ 



3 3 3 




u\s)u^{s) h' 



(s) +... 



i=0 j=0 fe=0 



3 



dh'js) 
ds 



g{A,K)=g{A,u) + eJ2M<s)) 



+ 




i=0 k=0 



When substituting these expansions into (2.1) we should remem- 
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ber that u is unit vector. Then for Sdei we get 



dh''{s) 



ds 



ds ■ 



Y 3 

S'def = S'-£ ^(mcuk{s) + - Ak{r{s)j^ 
i fc=o 



Let's apply integration by parts to first integral above. As a result 
we get the expression without derivatives of functions /i'^(s): 



fe=0 

Ve / (mcufe(s) + -^fc(r(s))) h^{s)ds- 



+ 



S2 



. fe=0 \ i=0 i=0 j=0 / 



Remember that function h^(^s) vanish at the ends of integration 
path h^{s\) = /i'^(s2) = (see §1). This provides vanishing of 
non-integral terms in the above formula for ^def- 

Now in order to derive differential equations for world line of 
particle we apply extremity condition for S. It means that term 
linear with respect to e in power expansion for 5def should be 
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identically zero irrespective to the choice of functions h''{s): 

1=0 1=0 j=0 

Let's calculate derivative in left hand side of (2.4). Then let's 
rearrange terms so that those with mc factor are in left hand 
side, while others with q/c factor are in right hand side: 




Now in right hand side of this equation wc find tensor of electro- 
magnetic field (see formula (9.9) in Chapter III). For transforming 
left hand side of this equation we use formula (11.3) from Chap- 
ter III. As a result we get the following equation for world line: 

\ i=0 i=0 / i=0 

In left hand side of the equation (2.5) we find covariant derivative 
with respect to parameter s along world line: 

3 

(2.6) mcWsUk=-y^Fkiu\ 

i=0 



Comparing (2.6) with the equations (11.9) from Chapter III, we 
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get formula for the vector of four-dimensional force acting on a 
particle with charge q in electromagnetic field: 

3 

(2.7) Ffe = ^^Ffe,M\ 

Suppose that we have rectangular Cartesian coordinate system 

in Minkowsky space. Then we can subdivide F into spatial and 
temporal parts and can calculate components of three-dimensional 
force vector: /' = ^J<? — \\^ (see formula (7.6) in Chap- 
ter III). Upon easy calculations with the use of formulas (7.2) and 
(9.4) from Chapter III for force vector f we get 

(2.8) f = 9E+^[v, H]. 

This formula (2.8) is exactly the same as formula for Lorentz force 
(see (4.4) in Chapter I). Thus formula (2.7) is four-dimensional 
generalization of formula for Lorentz force. Orthogonality condi- 
tion for 4-force and 4-velocity (see (7.5) in Chapter III) for (2.7) is 
fulfilled due to skew symmetry of tensor of electromagnetic field. 

Exercise 2.1. Prove that ga uge transformation of action func- 
tional (1.13) does not change dynamic equations of material point 
in electromagnetic field (2.6). 

Exercise 2.2. Verify that the relationship (7.5) from Chap- 
ter III holds for Lorentz force. 

§ 3. Dynamics of dust matter. 

Differential equation (2.6) describes motion of charged particles 
in electromagnetic field. If the number of particles is not large, 
then we can follow after the motion of each of them. When 
describing extremely large number of particles continual limit is 
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used, particles are replaced by continuous medium modeling their 
collective behavior. Simplest model describing large number of 
non-colliding particles is a model of dust cloud. In this model 




Fig. 3.1 Fig. 3.2 

particles of cloud move regularly (not chaotically). Their world 
lines can be modeled by regular family of lines filling the whole 
space (see Fig. 3.1). 

Another model is a model of ideal gas. Here particles also 
do not collide each other, i. e. their world lines do not intersect. 
However, their motion is chaotic (see Fig. 3.2). Therefore if we 
fill the whole space with their world lines, they would intersect. 

Besides two models considered mentioned above, there are 
models describing liquids and solid materials. Points of liquid 
and solid media move regularly (as on Fig. 3.1). However, in 
these media interaction of particles is essential. Therefore when 
describing such media one should either use detailed microscopic 
analysis and get macroscopic parameters by statistical averaging, 
or should use some heuristic assumptions based on experiment. 

In this book we consider only most simple model of dust cloud. 
In this case one should assume Minkowsky space to be filled by 
regular family of world lines. Some of them are world lines of real 
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particles, others are imaginary ones obtained by extrapolation in 
continual limit. Therefore at each point of M we have unit vector 
u, this is tangent vector to world line passing through this point. 
This means that dynamics of dust cloud can be described by 
vector field u(r). 

Apart from vector field u, below we need scalar parameter i/(r) 
which means the density of dust cloud. Wo define it as follows. 
Let's choose some small fragment of three-dimensional hypersur- 
face in M orthogonal to vector u(r) at the point r. The number 
of dust particles whose world linos cross this fragment is propor- 
tional to its three-dimensional volume: N = i/(r) V, parameter 
i/(r) is coefficient of proportionality. Parameter z^(r) has measure 
unit of concentration, it can be treated as concentration of par- 
ticles in small fragment of dust cloud near the point r measured 
in that inertial coordinate system for which particles of this small 
fragment arc at rest. By means of z^(r) and u(r) we compose new 
four-dimensional vector 



Vector (3.1) is called four-dimensional flow density for particles 
in cloud. If we choose some inertial coordinate system, then 
r]^/c is interpreted as concentration of particles in dust cloud, 
while other three components of four- dimensional vector 77 form 
three-dimensional vector of flow density. 

Suppose that dust cloud is formed by identical particles with 
mass m and charge q. Then four- dimensional current density 
vector can be represented as follows: 



By analogy with (3.2) one can define mass flow density vector. 



(3.1) 



r]{r) = ci^(r) u(r). 



(3.2) 



j(r) = qri{r). 



(3.3) 



/x(r) =m77(r). 
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Total number of particles in cloud is fixed. This conservation law 
is written as the following equality for rj: 

3 

(3.4) ^Vpr,P=0. 

p=0 

From (3.4) and (3.2) one can derive charge conservation law in 
form of the relationship (11.7) from Chapter III. Taking into 
account (3.3), we get rest mass conservation law: 

3 

(3.5) ^VpM^ = 0. 

p=0 

Rest mass conservation law here is fulfilled due to the absence of 
collisions when heavy particles can be produced from light ones 
(see § 7 in Chapter III). 

Let's consider dynamics of particles composing dust cloud. 
Vector field u is constituted by tangent vectors to world lines 
of dust particles. Therefore these world lines can be determined 
as integral curves of vector field u, i.e. by solving the following 
system of ordinary differential equations: 

(3.6) — =u'(r(s)), i = 0, ... ,3. 

Having determined world line of particle from differential equa- 
tions, we know vector of its 4-velocity u(s). Now let's calculate 
covariant derivative of vector u(s) with respect to parameter s: 

(3.7) w=^+j:i: u-{s). 

fe=0 71=0 
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Calculating derivative du^/ds in (3.7) we take into account (3.6) 
and the equality u(s) = u(r(s)). As a result we get 

ib=0 

Substituting (3.8) into (3.7), we derive the following formula: 

3 

(3.9) WsuP = ^u''VkuP. 

fc=o 

Right hand side of (3.9) is covariant derivative of vector field u(r) 
along itself (see more details in [3]). Substituting (3.9) into the 
equations of the dynamics of material point, we get: 

(3.10) VuU=— . 

mc 

Here F = F(r, u) is some external force field acting on particles 

of dust matter. For example in the case of charged dust in 
electromagnetic field the equation (3.10) looks like 



3 3 

(3.11) Y,^^y,u, = ^Y.P,,u\ 

fe=0 fe=0 

In contrast to the equations (11.9) from Chapter III, which 
describe dynamics of separate particle, here (3.10) are partial 
differential equations with respect to components of vector field 
u(r). They describe dynamics of dust cloud in continual limit. 
The equation for scalar field u{r) is derived from conservation law 
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(3.4) for the number of particles. Combining these two equations, 
we obtain a system of differential equations: 

mc 

fe=o 

3 3 
fe=0 fc=0 

System of partial differential equations (3.12) yields complete 
description for the dynamics of dust cloud. 

Model of dust matter cam be generalized a little bit. We can 
consider mixture of particles of different sorts. For each sort of 
particles we define its own vector field u(i,r) and its own scalar 
field of concentration J^(i,r). Then formulas (3.2) and (3.3) for j 
and n are generalized as follows: 

n n 

= Y1 = '^(*) '^)- 

i=l i=l 

Here 'q{i, r) = cu{i, r) u(i, r). Each pair of fields u(i, r) and y(i, r) 
satisfies differential equations (3.12). We can derive mass and 
charge conservation laws from these equations. 

§ 4. Action functional for dust matter. 

Let's study the dynamics of dust matter in electromagnetic 
field within framework of Lagrangian formalism. Fort this purpose 
we need to pass to continual limit in action functional (1.8). 
For the sake of simplicity we consider dust cloud with identical 
particles. Omitting details of how it was derived, now we write 
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(3.12) 
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ultimate formula for action functional (1.8) in continual limit: 

V2 

S = -m j ^/g{r],r])^/-detg d^r- 

Vi 

V2 

(4-1) J g^rj,A)^/^d^d''r- 

Vi 

V2 3 3 

p=Ofc=0 

Instead of deriving formula (4.1) from (1.8) we shall verify this 
formula indirectly. For this purpose we shall derive dynamical 
equation (3.11) from principle of extremal action applied to action 
functional (4.1). 

For describing dust matter in (4.1) we have chosen vector field 
T7(r) defined in (3.1). Other two fields u(r) and p{r) can be 
expressed though vector field T]{r): 



(4.2) cu = I77I = y/gi'n,'n) , u= — . 

c 1^ 

Dealing with variation of vector field 17 (r) we should always 
remember that components of this field are not independent 
functions. They satisfy differential equation (3.4). In order 
to resolve tis equation (3.4) we use slightly modified version of 
theorem 10.1 from Chapter III. 

Theorem 4.1. Let M be some n-dimensional manifold, where 
n>2, equipped with metric gij . For each vector field rj with zero 
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divergency with respect to metric connection 

n 

(4.3) ^Vp7?f = 

p=i 

there is skew-symmetric tensor Geld (p of type (2, 0) such that the 
following relationships are fulfilled 

n 

(4.4) ^ = ^V,v^«. 

9=1 

Proof. Writing relationships (4.3), we use well-known for- 
mula for components of metric connection, see formula (11.3) in 
Chapter III. As a result we get 

P=l p=l p=l 5=1 p=l 

( dcjpk dgks dgps \ ^ 
^ 21^1^1^^ \ dr" drP dr'' J ' ' 

p=l s=l fc=l ^ ^ 

Note that last two derivatives of metric tensor in round brackets 
are canceled when we sum over indices p and k. This is because 
gP'^ is symmetric. Hence 

E%.^ = EgF + ^EEE^^'^|?^^ = 

P=l p=l p=l 5=1 fc=l 

p=i p=i ^ ' 
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For further transforming of this expression (4.5) we use well 
known formula for logarithmic derivative of determinant: 

, , 9 In I det q\ f t do 

Substituting (4.6) into (4.5), we transform (4.5) so that 



(47) 1 y- gfa^Vpe^) 

Let's carry out analogous calculations for right hand side of (4.4) 
taking into account skew symmetry of the field ip^'' and symmetry 
of connection components T^^. These calculations yield 



(48, ^v,^.. ' ^^'^"M. 



Denote i'P = ^ I dot g \ rp and V-"^' = \/ \ dot g \ ^P'^ . Now on the 
base of (4.7) and (4.8) it is easy to understand that proof of 
theorem 4.1 is reduced to theorem 10.1 from Chapter III. □ 

Remark. Generally speaking, theorem 10.2 has no direct 

generalization for the case of spaces with metric. It is generalized 
only for metric spaces with zero curvature tensor R\.^^ = 0. 

Let's define deformation of the field r; in a way similar to that 
we used for vector-potential A in § 1: 

(4.9) r}f(r)=r7f(r)+£C^(r) + ... . 



Both fields 17 and r] satisfy differential equation (3.4). Hence 
vector field C, defined in (4.9) also satisiy this equation. Let's 
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apply theorem 4.1 to vector field 

3 

(4.10) C^ = ^Vfe<^^"=. 

fc=0 

Theorem 4.1 does not specify tensor field ifp'^ in (4.10), this can 
be any skew-symmetric tensor field. However, we choose it in very 
special form as follows: 

(4.11) ifiP^ = rf h!" - hF T]^ . 

This choice can be motivated by the following theorem. 

Theorem 4.2. For any two vector fields C and t], where rj 0, 
both satisfying differential equation (3.4) tiiere is vector field h 
such that vector field ^ is given by formula 

3 

C^ = 5^Vfc(r/'/i'=-/ifj7'=). 

fc=0 

Our choice (4.11) leads to the following expression for the field r): 

3 

(4.12) f]P{r) = r,P{r) + e ^ Vkiv" h'' - hP v"") + ■ ■ ■ ■ 

fe=o 

Quantities /i*(r) in (4.12) are chosen to be smooth functions being 
nonzero only within some restricted domain CI in Minkowsky 
space. 

When substituting (4.12) into action functional (4.1) we use 
the following expansion for \/g{i), 17) : 
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We have analogous power expansion for the expression under 
second integral in formula (4.1): 

3 3 

g(»7,A) =5(J7, A) +e^^^pVfc^f '= + ... . 

p=0 fc=0 

Substituting these two expansions into (4.1), we take into account 
(4.2). For the action 5def this yields 

3 3 



Y^Y^Up Vfc<pfV-det5rfV- 

O P=Ok=0 

/3 3 
^ ^ Vfc¥'P'= y^d^rfV + . . . . 
^—n I — n 



(4-13) 

eq 



p=0 k=0 

Further in order to transform the above expression (4.13) we 
use Ostrogradsky-Gauss formula. In the space equipped with 
metric this formula is written as follows: 

3 

(4.14) l'Y,ykz''V^d^d\= f g{z,n)dV. 

fe=o 



an 



Here z^, z^, z^ are components of smooth vector field z, 
while n is unit normal vector for the boundary of the domain 
n. In order to transform first integral in formula (4.13) we take 
z'' = J2p^Q Up ipP^ . Then in right hand side of (4.14) we obtain 

3 3 3 3 3 

/s=0 p=0 fe=0 p=0 k=0 
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Right hand side of (4.14) vanishes since </5^''^ do vanish on the 
boundary of ft. Hence we have the equality 



/3 3 
a P=Ofc=o 

3 3 
J ^ — n ! — n 



In a similar way we transform second integral in (4.13). In whole 
for the action ^def we get the following expression 



(4.15) 



Sdef = S + em I Y^Y^ VkUp ifP'' ^/ - det g d^r + 

p=0 fe=0 



3 3 



J n l._r, 



p=0 /c=0 



Extremity of action S means that linear part with respect to e 
in formula (4.15) should vanish: 



3 3 



Y Y {^^kUp + ^VkAp) v'fV-det^dV = 0. 

Let's substitute formula (4.11) for t/j^''^ into the above equality. 
Then it is transformed to the following one 

/3 3 
Y Y {^^kUp + ^VkAp) r/P/i'= ^-deig d\ = 
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/3 3 

Let's exchange indices k and p in second integral. Thereafter 
integrals can be united into one integral: 

/3 3 

(4.16) n fe=Op=o ^ 

-^^pAk) ?7^/iV-det5dV = 0. 

Now let's take into account that in resulting equality = h''{r) 

are arbitrary smooth functions vanishing outside the domain Q. 
Therefore vanishing of integral (4.16) means vanishing of each 
summand in sum over index k in the expression under integration: 

3 

(4.17) (m WkUp - m VpUk + \ i^fcp) rf = Q. 

Here we used the relationship (11.5) from Chapter III. It re- 
lates tensor of electromagnetic field and four-dimensional vector- 
potential. 

In order to bring the equation (4.17) just derived to its ultimate 
form we use the relationships (4.2), which relate vector field 77 
and vector field u: rf = cuu^. Since u is unit vector, we have 

3 

(4.18) ^ufVfcUp = 0. 

p=0 
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Taking into account (4.18), we bring (4.17) to the following form: 



Now it is easy to see that (4.19) exactly coincides with the equa- 
tion (3.11), which we have derived earlier. This result approves 
the use of the action (4.1) for describing charged dust matter in 
electromagnetic field. 

Exercise 4.1. Prove that for any skew-symmetric tensor field 
ifP"^ vector field rj determined by formula (4.4) has zero divergency, 
i. e. is satisfies differential equation (3.4). 

Exercise 4.2. Prove theorem 4.2. For this purpose use the 
following fact known as theorem on rectification of vector field. 

Theorem 4.3. For any vector Geld rj ^ there exists some 
curvilinear coordinate system r°, r^, r^, such that r]^ = 1, 
T]^ = 0, r]"^ = 0, 7]^ = in this coordinate system. 

Exercise 4.3. Prove theorem 4.3 on rectification of vector field. 

Exercise 4.4. Derive Ostrogradsky-Gauss formula (4.14) for 
the space equipped with metric on the base of the following integral 
relationship in standard space R": 



3 



3 



(4.19) 





n 



on 



§ 5. Equations for electromagnetic field. 



In this section we continue studying action functional (4.1). 
This functional describes dust cloud composed of particles with 
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mass m and charge q in electromagnetic field. In previous section 
we have found that applying extremal action principle to S with 
respect to the field 77 one can derive dynamical equations for 
velocity field in dust cloud. Now we shall apply extremal action 
principle to S with respect to vector-potential A. Deformation of 
vector-potential is defined according to (1.9), (1.10), (1-11), (1-12): 

(5.1) Ai{r) = Ai{r) + ehi{r) + . . . . 

For components of tensor of electromagnetic field we derive 

(5.2) Fij = Fij + e {S/ihj - Vjhi) + ... . 

When substituting (5.2) into action functional (4.1) we carry out 
the following calculations: 

tt^^'^^'' = tttt hk 4- 5- = 

p=0 fe=0 i=0 j=0 p=0 fe=0 

3 3 3 3 

= E E ^pfc + 2e^ ^Ff'=(Vp/ife - VkK) + .... 

p=0 k=0 p=0 fc=0 

Taking into account skew symmetry of tensor F^'', this expansion 
can be simplified more and can be brought to the form 

33 33 33 

^ ^ Fpfe = ^ ^ Fpfe + ^ Ff*=Vj,/ifc + . . . . 

p=0 fe=0 p=0 fe=0 p=0 fe=0 

Analogous calculations in substituting (5.1) into (4.1) yield 

3 

g{r],A)=g{r],A)+e^ri''hk + ... ■ 

k=0 
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As a result for deformation of action functional (4.1) we get 

. 3 



n 



fe=0 



£ 

Attc 



/3 3 
^ — n I — n 



p=0 fe=0 



Let's transform second integral in the above expansion for 5def 
by means of Ostrogradsky-Gauss formula (4.14). For this purpose 
let's choose = J2k=o ^^^^k and take into account vanishing of 
hk on the boundary of the domain O. Then for 5def we get 

Q fe=0 \ p=0 / 

Extremal action principle means that linear in e part of the above 
expansion for S'dof should vanish. Note also that is an arbitrary 
domain and /ife(r) are arbitrary functions within fl. This yield 
the following equations for the tensor of electromagnetic field: 

(5.3) ^VpFPfe = l^^'=. 

p=0 

Remember that ri{r) is related to current density by means of 
(3.2). Then (5.3) can be written as 

(5.4) ^VpF'=f = -^/. 

p=0 
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It is easy to see that (5.4) are exactly Maxwell equations written 
in four-dimensional form (see (11.4) in Chapter III). Another pair 
of Maxwell equations written in four-dimensional form 

g=0 fe=0 s=0 

is a consequence of the relationship Fpq = VpAq — VqAp (see 
formula (11.5) in Chapter III). 

Exercise 5.1. Which form will have differential equations (5.3) 
if we consider dust cloud composed by particles of several sorts 
with muNscs rn{l), . . . , m{N) and charges q{l), . . . , q{N) ? Will 
this change differential equations (5.4) ? 



CHAPTER V 



GENERAL THEORY OF RELATIVITY 



§ 1. Transition to non-flat metrics 
and curved Minkowsky space. 

Passing from classical electrodynamics to special theory of re- 
lativity, in previous two chapters we have successively geometrized 
many basic physical concepts. Having denoted r" = ct and com- 
bining with components of three-dimensional radius-vector in 
inertial coordinate system, we have constructed four-dimcnsional 
space of events (Minkowsky space). This space appears to be 
equipped with metric of signature (1,3), which is called Min- 
kowsky metric. Thereby inertial coordinate systems are inter- 
preted as orthonormal bases in Minkowsky metric. 

In four-dimensional formalism dynamics of material point is 
described by vectorial differential equations, while Maxwell equa- 
tions for electromagnetic field are written in tensorial form. Due 
to this circumstance in previous two chapters we managed to 
include into consideration skew-angular and even curvilinear co- 
ordinate systems in Minkowsky space. Thereby we got explicit 
entries of metric tensor components Qij, metric connection com- 
ponents , and covariant derivatives in all our equations. 

Next step in this direction is quite natural. One should keep 
the shape of all equations and pass from flat Minkowsky metric 
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to metric of signature (1, 3) with nonzero curvature tensor: 

npfe npA- 3 3 

(-[ -i\ Ti^ = J'? _ ''I J. \ " pfe ps _ \ " pfe ps 

s=0 s=0 

This crucial step was first made by Einstein. Theory he had dis- 
covered in this way later was called Einstein's theory of gravitation 
or general theory of relativity. 

Definition 1.1. Four- dimensional afRne space equipped with 
orientation, polarization, and with metric of signature (1, 3) and 
nonzero curvature (1.1) is called curved Minkowsky space. 

In non-flat Minkowsky space we loose some structures available 
in flat case. In such space there are no coordinates for which 
Minkowsky metric is given by matrix (2.7) from Chapter III, i. e. 
here wc have no incrtial coordinate systems. This is substantial 
loss, but it is not catastrophic since dynamic equation for material 
points and Maxwell equations rewritten in vectorial and tensorial 
form arc not bound to inertial coordinate systems. 

Geodesic lines in curved Minkowsky space do not coincide 
with afiine straight lines. Therefore afiine structure becomes 
excessive restriction in general relativity. As appears, one can 
give up topologic structure of flat space M** as well. Even in 
two-dimensional case, as we know, apart from deformed (curved) 
plain, there arc surfaces with more complicated topology: sphere, 
torus and sphere with several handles glued to it (see [5]). In 
multidimensional case these objects are generalized in concept of 
smooth manifold (see details in [2], [5], and [6]). 

Smooth manifold M of dimension n is a topologic space each 
point of which has a neighborhood (a chart) identical to some 
neighborhood of a point in R". In other words M is covered by a 
family of charts Ua, each of which is diffeomorphic to some open 
set Va in K". Such chart maps (chart diffeomorphisms) define 
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local curvilinear coordinate systems within their chart domains 
Ua- At those points of manifold M, where two chart domains are 
overlapping, transition functions arise. They relate one curvilinear 
coordinate system with another: 



f' = r'(r^, . . . , r"), where i = l,.. . , n, 
r' = r'(f^, . . . , f"), where i = 1, . . . , n. 



According to definition of smooth manifold, transition functions 
(1.2) are smooth functions (of class C°°). Transition functions 
determine transition matrices S and T: 



Presence of transition matrices (1.3) lead to full-scale theory of 
tensors, which is almost literally the same as theory of tensors 
for curvilinear coordinates in R" (see [3]). The only difference 
is that here we cannot choose Cartesian coordinates at all. This 
is because in general there is no smooth diffeomorphic map from 
manifold M to M". 

Definition 1.1. Four-dimensional smooth manifold equipped 

with orientation, polarization, and with metric of signature (1,3) 
is called generalized Minkowsky space or Minkowsky manifold. 

§ 2. Action for gravitational field. 
Einstein equation. 

Space of events in general relativity is some smooth Minkowsky 
manifold M. This circumstance provides additional arbitrariness 
consisting in choosing M and in choosing metric on M . Nonzero 
curvature described by tensor (1.1) is interpreted as gravitational 
field. Gravitational field acts upon material bodies and upon 
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electromagnetic field enclosed within M. This action is due to 
the presence of covariant derivatives in dynamic equations. The 
magnitude of gravitational field itself should be determined by 
presence of matter in M in form of massive particles or in form of 
electromagnetic radiation, i. e. we should have backward relation 
between geometry of the space and its content. 

In order to describe backward relation between gravitational 
field and other physical fields we use Lagrangian formalism and 
extremal action principle. Let's start from action functional (4.1) 
in Chapter IV. It is sum of three integrals: 

(2.1) S = Smat + Sint + Se\. 

First integral ^mat is responsible for material particles in form of 
dust cloud, second integral describes interaction of dust cloud and 
electromagnetic field, third term in (2.1) describes electromagnetic 
field itself. In order to describe gravitational field one more 
summand in (2.1) is added: 

(2.2) S = Sgi + Smat + Sint + Se\. 

This additional term is chosen in the following form: 



Here 7 gravitational constant same as in Newton's universal law 
of gravitation (see formula (1.11) in Chapter I). Scalar quantity R 
in (2.3) is scalar curvature determined by curvature tensor (1.1) 
according to the following formula: 



(2.3) 




3 



3 



3 



(2.4) 




q=0 fe=0 j=0 
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Ricci tensor is an intermediate object relating curvature tensor 
(1.1) and scalar quantity (2.4). Here are its components: 



Ricci tensor is symmetric (see [3]). Scalar curvature R is obtained 

by contracting Ricci tensor and metric tensor g'^^ with respect to 
both indices q and j. This fact is obvious due to (2.5) and (2.4). 

Note that sometimes in the action for gravitational field (2.3) 
one more constant parameter A is added: 



This parameter is called cosmological constant. However, accord- 
ing to contemporary experimental data the value of this constant 
is undetectably small or maybe is exactly equal to zero. Therefore 
further we shall use action in form of (2.3). 

Note also that metric tensor describing gravitational field enters 
in implicit form into all summand in (2.2). Therefore we need 
not add special terms describing interaction of gravitational field 
with material particles and electromagnetic field. Moreover, such 
additional terms could change the form of dynamical equations for 
matter and form of Maxwell equations for electromagnetic field 
thus contradicting our claim that these equations are the same in 
general and in special relativity. 

Now let's begin with deriving dynamical equations for grav- 
itational field. For this purpose we consider deformation of 
components of metric tensor given by the following relationship: 



3 



(2.5) 




fe=0 




{R +2 A) ^y - det g d^r. 



(2.6) 



g'\r)=g'^{r)+eh'^{r) + . 
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Functions h''^{r) in (2.6) are assumed to be smooth functions 
vanishing outside some restricted domain Q c M. Deformation of 
matrix g^^ lead to deformation of inverse matrix gij : 

9ij ~ 9ij ^ ^ij ~^ ' • • — 

(2.7) 3 3 

= 9ij - ^Yli^Siph^'^ 9qj + ■ ■ ■ . 

Let's differentiate the relationship (2.7) and let's express partial 
derivatives through covariant derivatives in resulting formula: 

Qj.k Qj-k Qj-k Qj-k '-J 

(2.8) 3 3 

In (2.8) wo used covariant derivatives corresponding to non- 
deformed metric gij. Now on the base of (2.8) we calculate 
the following combination of derivatives: 

dgkj , dcjik dcjij _ dgkj dgik dgij 



Qfi Qj-j Qj-k Qj.t Qj.j Qj.k 

(2.9) / 3 X 

- e I Vj/ifcj + Vjhik - Vkhij -2^T^jhpk\ + ... . 

\ p=0 / 

Let's use the relationships (2.6) and (2.8) in calculating defor- 
mation of connection components. For this purpose let's apply 
well-known formula to T^j (see formula (11.3) in Chapter III): 

3 
fe=0 
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This expansion for T^j can be written in symbolic concise form 

(2.10) t^.=^+eYi; + ... 

by introducing the following quite natural notation: 

1 3 

(2.11) = 2 E a""" (^i^kj + yjhik - Vkhij) . 

fe=0 

Now let's substitute the expansion (2.10) into the formula (1.1) 
for curvature tensor. This yields 

(2.12) = + e (V,l}^, - V.-F^) + . . . . 

Upon contracting (2.12) with respect to one pair of indices we get 
similar expansion for deformation of Ricci tensor: 

3 

(2.13) R,^ = R,j + e J2 C^kYjl - V.F^) + . . . . 

fe=0 

Wc multiply (2.13) by gqj using formula (2.6). Then wc carry out 
complete contraction with respect to both indices q and j. This 
yields deformation of scalar curvature: 

3 3/ 3 \ 

^ = + £ ^ ii,,- /i« + ^ ( Vfc Y]l - V,- f4 ) + . . . . 

j=0 9=0 \ fe=0 / 

Let's introduce vector field with the following components: 

j=0 g=0 
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Then we can rewrite deformation of scalar curvature R as 

3 3 3 

(2.14) i? = i? + e^^i?gj/i«+£^Vfe^'= + ... . 

j=0 q=0 fc=0 

When substituting (2.14) into action integral (2.3) wc should note 
that second sum in (2.14) is exactly covariant divergency of vector 
field Z. Components of Z are smooth functions vanishing outside 
the domain Cl. Therefore integral of such sum is equal to zero: 

3 

j J2 ^kZ^ V- det g d'^T = j g{Z, n) dV = 0. 
n '=='^ an 

This follows from Ostrogradsky-Gauss formula (see (4.14) in 
Chapter IV). Hence for deformation of we get 



In deriving this formula we also used the following expansion: 

(2.15) v/^d^ = v/^d^ f 1 - ^ E E + • • • • 

It follows from (2.6). Now we shall not calculate deformations of 
other three terms in (2.2) in explicit form. This will be done in 
§ 4 and § 5 below. However, we introduce notation 

(2.16) 

Here S'm.f. denotes overall action for all material fields other than 
gravitation. The number of terms in the sum (2.16) could be much 
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more than three, if one consider more complicated models for 
describing matter. But in any case action of gravitational field is 
excluded from this sum since gravitational field plays exceptional 
role in general relativity. Now we shall write deformation of the 
action (2.16) in the following conditional form: 



(2.17) 5def = S^.i. + y Y.Y.^oi V-detg d^r + . . . . 



Then extremity condition for total action (2.2) is written as 



This equation (2.18) is known as Einstein equation. It is ba- 
sic equation describing dynamics of metric tensor gij in general 
theory of relativity. 

Exercise 2.1. Derive the relationships (2.7) and (2.15) from 
the expansion (2.6) for deformation of tensor g^^ . 



Tensor T in right hand side of Einstein equation (2.18) is called 
energy-momentum tensor for material fields. It is determined 
by the relationship (2.17) and comprises contributions from all 
material fields and their interactions. In the model of dust matter 
in electromagnetic field tensor T is composed of three parts (see 
formula (2.16)). 

Energy-momentum tensor is related with 4-momentum conser- 
vation law for material fields. In order to derive this conservation 
law we use well-known Bianchi identity: 




3 



3 



(2.18) 




§ 3. Four-dimensional momentum 
conservation law for fields. 



(3.1) 



154 



CHAPTER V. GENERAL RELATIVITY. 



More details concerning Bianchi identity (3.1) can be found in [2] 
and [6]. Let's contract this identity with respect to i and p: 



3 



(3.2) VkRsj + 2^ Vpi?^., - VjR.k = 0. 

p=0 

Here we used skew symmetry of curvature tensor with respect 
to last pair of indices (see [3]). Let's multiply (3.2) by g^^ and 
contract it with respect to double indices s and j. Upon slight 



transformation based on skew symmetry R^J = —R^j we get 

3 



(3.3) ^ VsRi - iVfci? = 0. 

Now let's raise index j in the equation (2.18), then apply covariant 
differentiation Vj and contract with respect to double index j: 

(3-4) ^:v,^^^v,i^=^x:v.•^•• 

j=0 j=0 

Comparing (3.3) and (3.4), we get the following equation for 
energy-momentum tensor of material fields: 

3 

(3.5) ^V,T,^=0. 

3=0 

The equation (3.5) expresses A-momentum conservation law for 
the whole variety of material fields. It is usually written in the 
following form with raised index q: 

3 

(3.6) $^V,T«=0. 

Energy-momentum tensor is symmetric therefore the order of 
indices q and j in (3.6) is unessential. 
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§ 4. Energy-momentum tensor 
for electromagnetic field. 

Energy-momentum tensor for whole variety of material fields 
is defined by the relationship (2.17). By analogy with (2.17) we 
define energy-momentum tensor for electromagnetic field: 

/3 3 
^ J2 T,j /i«v^d^ rfV + . . . . 
n 3=0 

Basic fields in the action 5ei are covariant components of vector- 
potential Ai{r). Covariant components of tensor of electromag- 
netic field are defined by formula 

(4.2) F,, = V.A, - y,A, = ^-^ 

(sec also formula (11.5) in Chapter III). Ultimate expression in 
right hand side of (4.2) has no entry of connection components 
r^^. Therefore covariant components Fij are not changed by 
deformation of metric (2.6). Upon raising indices we get 

3 3 

and, using this formula, for contravariant components F^i of 
tensor of electromagnetic field we derive the expansion 

3 3 

(4.3) FP*^ = + £ ^ ^{h^' 9^^ + 9^" h!"^ ) Fij + . . . . 

j=0 j=0 
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Substituting F^'^ and g into action functional Se\, we get 

V2 



3 3 



p=Ofe=0 

Then, taking into account (4.3) and (2.15), we derive formula 
5aef = 5e, - ^ / ^ E f E E 2 F,, F,, - 

^2 q=0 i=0 \ p=0 i=0 

-2EE^f^^'"5« P''V^d^'^'^ + --- • 

p=0 i=0 / 

Comparing this actual expansion with expected expansion (4.1) 
for Sdei, we find components of energy-momentum tensor for 
electromagnetic field in explicit form: 

(4-4) T,, = -^EE (Ppia^'Fij - \Fp^F'" 



3 3 
p=0 i=0 

Raising indices q and j in (4.4), for contravariant components of 
energy-momentum tensor T we derive 

33, 

(4.5) = - ^ E E ( - i^p^ 

p=0 i=0 ^ 

By means of formula (4.5) one can calculate covariant divergency 
for energy-momentum tensor of electromagnetic field: 
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Formula (4.6) shows that 4-momentum conservation law for sepa- 
rate electromagnetic field is not fulfilled. This is due to momen- 
tum exchange between electromagnetic field and other forms of 
matter, e. g. dust matter. 

Exercise 4.1. Verify the relationship (4.6). For this purpose 
use well-known formula for commutator of covariant derivatives 

3 

{ViVJ-V^V,)Ak = -J2^kiJAs 
s=0 

and properties of curvature tensor (see details in [3]). 

Exercise 4.2. Calculate components of energy-momentum 
tensor (4.5) in inertial coordinate system for Eat Minkowsky met- 
ric. Compare them with components of Maxwell tensor, with den- 
sity of energy, and with vector of energy Bow for electromagnetic 
field (see formulas (2.5) and (2.15) in Chapter II). 

§ 5. Energy-momentum tensor 
for dust matter. 

Let's consider energy- momentum tensor related with last two 
terms Smat and ^int in the action (2.16). They contain entries of 
vector field rf whose components satisfy differential equation 

3 

(5.1) ^Vpr,^=0, 

see (3.4) in Chapter IV. This circumstance differs them from com- 
ponents of vector-potential A. Metric tensor gij enters diff'erential 
equation (5.1) through connection components T^j of metric con- 
nection. Therefore by deformation of metric gij gij one cannot 
treat ry^ as metric independent quantities. 
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In order to find truly metric independent variables for dust 
matter we use formula (4.7) from Chapter IV and rewrite differ- 
ential equation (5.1) as follows: 



Denote ff ~ "rf \J — dot g . These quantities ff can be treated as 
metric independent ones since differential constraint for them is 
written in form of the equation that does not contain metric: 



Expressing rf through for action functional ^int describing 
interaction of dust matter and electromagnetic field we get 



It is easy to see that integral (5.3) does not depend on metric 
tensor. Therefore action functional S'int makes no contribution to 
overall energy-momentum tensor. 

Now let's express ry^ through if in action functional Smat for 
dust matter. As a result we get formula 



p=0 




(5.2) 




p=0 



(5.3) 




(5.4) 



Smut = -m 




The dependence of this functional on metric tensor is completely 
determined by explicit entry of gpq under square root sign in right 
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hand side of (5.4). Therefore power extension for Sma,t is easily 
calculated on the base of the expansion (2.7): 

Let's compare this expansion with expected expansion for 5def: 

.33 

5def = 5mat + f / E E /^"'v^-detfl + . . . . 

n n 



n 



p=0 g=0 



By this comparison we find explicit formula for components of 
energy-momentum tensor for dust matter: 



(5.5) Tpq = ''P 'li = jjic \/g{r],r]) Up Uq. 

Contravariant components of energy-momentum tensor (5.5) 
are obtained by raising indices p and q: 

mc rf rf 



(5.6) Tf« = = ™c ^/g^ 

Using collinearity of vectors u and r/ (sec formula (3.1) in Chap- 
ter IV) and recalling that u is unit vector, we can bring formula 

(5.6) to the following simple form: 

(5.7) Tf*^ = rncwf^7^ 

Formula (5.7) is convenient for calculating covariant divergency of 
energy-momentum tensor for dust matter: 

3 3 

(5.8) E'^^^'" = f E^"'''- 

s=0 s=0 
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Now, applying formula (3.2) from Chapter IV, we can transform 
formula (5.8) and write it as follows: 



Let's compare (5.9) with analogous formula (4.6) for energy- 
momentum tensor of electromagnetic field. Right hand sides of 
these two formulas differ only in sign. This fact has transparent 
interpretation. It means that in our model the overall energy- 
momentum tensor for matter 



satisfies differential equation (3.6). This fact is in complete 
agreement with 4-momentum conservation law. 

Another important conclusion, which follows from of (4.6) and 
(5.9), is that 4-momentum conservation law for the whole variety 
of material fields can be derived from dynamical equations for 
these fields. Therefore this law is valid also in special relativity, 
where Einstein equation (2.18) is not considered and where in 
general case for fiat Minkowsky metric it is not fulfilled. 

Exercise 5.1. Derive the relationship (5.8) on the base of equa- 
tions (3.4) and (4.19) from Chapter IV. 



Event space in general theory of relativity is some Minkowsky 
manifold M with Minkowsky metric of signature (1,3). This 
metric is determined by material content of the space according 
to Einstein equation (2.18). However, topology of the manifold 
M has great deal of arbitrariness. This manifold can have local 
singularities at the points with extremely high concentration of 



(5.9) 




Tm.f. — Tmat + Tel 
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matter. Such objects are called black holes. Moreover, global 
topology of M also can be nontrivial (other than topology of R^). 
In contemporary physics most popular models of M include big 
bang in the very beginning of times. According to these models 
in far past times our Universe M was extremely small, while 
density of matter in it was extremely high. In further evolution 
our Universe was expanding up to its prc!sciit size. Will this 
expansion last infinitely long or it will change for contraction? 
This problem is not yet solved. The answer to this question 
depends on estimates of total amount of matter in the Universe. 

In this book wc cannot consider all these fascinating problems 
of modern astrophysics and cosmology. However, I think the 
above theoretical material makes sufficient background for to 
continue studying these problems e.g. in books [2], [7], and [8]. 
I would like also to recommend the book [9] of popular genre, 
where these problems are discussed in commonly understandable 
and intriguing manner. 
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